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Abstract. We consider two models for directed polymers in space-time independent ran- 
dom media (the O'Connell-Yor semi-discrete directed polymer and the continuum directed 
random polymer) at positive temperature and prove their KPZ universality via asymptotic 
analysis of exact Fredholm determinant formulas for the Laplace transform of their partition 
functions. In particular, we show that for large time r, the probability distributions for 
the free energy fluctuations, when rescaled by r^/'^, converges to the GUE Tracy- Widom 
distribution. 

We also consider the effect of boundary perturbations to the quenched random media on 
the limiting free energy statistics. For the semi-discrete directed polymer, when the drifts of 
a finite number of the Brownian motions forming the quenched random media are critically 
tuned, the statistics are instead governed by the limiting Baik-Ben Arous-Peche distributions 
from spiked random matrix theory. For the continuum polymer, the boundary perturbations 
correspond to choosing the initial data for the stochastic heat equation from a particular class, 
and likewise for its logarithm ~ the Kardar-Parisi-Zhang equation. The Laplace transform 
formula we prove can be inverted to give the one-point probability distribution of the solution 
to these stochastic PDEs for the class of initial data. 
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1. Introduction and main results 

The main results of this paper are contained in Sections II .21 II .31 and II .4[ The below section 
introduces the study of directed polymers and motivates our main results within that field. 

1.1. Directed polymers in random media. This article studies the effect of quenched 
disorder (homogeneous and inhomogeneous) on a class of path measures introduced first 
by Huse and Henley [37] which are commonly called directed polymers in a random media 
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(DPRM). Such polymers are directed in what is often referred to as a time direction, and 
then are free to configure themselves in the remaining d spatial dimensions. The probability 
dPp^Q^TT^-)) of a given configuration 7r(-) of the polymer is then given relative to an underlying 
reference path measure dPo by a Radon-Nikodym derivative. This derivative is written as a 
Boltzmann weight involving a Hamiltonian Hq which assigns an energy to the pathQ: 



In the above equation (3 is the inverse temperature which balances the entropy of the under- 
lying reference path measure with the energy of the Hamiltonian. The subscript Q stands 
for quenched which means that this Hq{7t{-)) is a function of some disorder u which we think 
of as an element of a probability space. Hence, with respect to this probability space, Hq is 
a random function. The normalization constant Z^.q, given by 



is the quenched partition function, and is a function of uj as well. The measure dPp^q is a 
quenched polymer measure since it is also a function of u. We denote averages with respect 
to the disorder u by E, so that E(Z^ g) is the average of the quenched partition function. We 
denote by P the probability measure for the disorder uj and denote the variance with respect 
to the disorder as Var -. From now on we will assume dP^ is the path measure of a random 
walk with either a free end point, or a specified pinned endpoint. The latter case is called a 
point-to-point polymer. We will focus mainly on point-to-point polymers herein. 

At infinite temperature, /3 = 0, and under standard hypotheses on dPo (i.e., i.i.d. finite 
variance increments) the measure ciP/j^Q (7r(-)) rescales diffusively to that of a Brownian mo- 
tion and thus the polymer is purely maximizing entropy. At zero temperature, /3 = oo, the 
polymer measure concentrates on the path (or paths) vr which maximize the polymer energy 
Hq{t[). a well studied challenge is to understand the effect of quenched disorder at positive 
/3 on the behavior of a dP^^g-typical path of the free energy Fq^p := /3~^ln(Z^ Q). A rough 
description of the behaviour is given by the transversal fluctuation exponent and the lon- 
gitudinal fluctuation exponent x- There are many different ways these exponents have been 
defined, and it is not at all obvious that they exist for a typical polymer model - though it 
is believed that they do. As n goes to infinity, the first exponent describes the fluctuations 
of the endpoint of the path vr: typically |7r(?T.)| ^ n^. The second exponent describes the 
fluctuations of the free energy: VarF^^g ^ n^^. Assuming the existence of these exponents, 
in order to have a better understanding of the system it is of essential interest to understand 
the statistics for the properly scaled location of the endpoint and fluctuations of the free 
energy. 

From now on we will focus entirely on Hamiltonians which take the form of a path inte- 
gral through a space-time independent noise field (the quenched disorder). In the discrete 
setting of dPo given as a simple symmetric random walk (SSRW) of length n, we consider a 
homogeneous noise field given by i.i.d. random variables Wt^x- The quenched Hamiltonian is 
then defined as Hq{'k{-)) = Ylt=o^t,n{t)- 



Actually, the energy is —Hq{tt{-)) but we keep the same convention as in the models analyzed. 



= ^,"ie/^^«(-«)rfPo(vr(-)). 
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The first rigorous mathematical work on directed polymers was by Imbrie and Spencer [IQ] 
in 1988 where they proved that in dimensions d > 3 and with small enough /3, the polymer 
is diffusive = 1/2). Bolthausen [15] strengthened the result to a central limit theorem 
for the endpoint of the walk. This means that in c? > 3 entropy dominates at high enough 
temperature, i.e., the polymer behaves as if there were no noise. Bolthausen's work relied on 
the now fundamental observation that renormalized partition function (for (IPq a SSRW of 
length n) Wn = Zj^^Q/Ki^Zp^o) is a martingale. By a zero-one law, the limit W^o = lim„_j.oo 
is either almost surely or almost surely positive. 

Since at /3 = the effect of the randomness u vanishes and one has Wqo = 1, one refers 
to the case Woo > as weak disorder, while the term strong disorder is used when Woo = 0. 
There is a critical value Pc such that weak disorder holds for [i < Pc and strong for /3 > (5c- 
It is known that /3c = for d G {1, 2} [21] and < /3c < cxo for d > 2>. In (i > 3 and weak 
disorder the walk converges to a Brownian motion [22]. On the other hand, the behavior of 
the polymer paths in strong disorder is different since there exist (random) points at which 
the path vr has a uniform (in n) positive probability (under rfPg q) of ending (see [2U|[2T]). 

The KPZ universality conjecture for d = 1. In this paper we focus on the d = 1 case. The 
universality conjecture says that the scaling exponents and limiting fluctuation statistics for 
the free energy exist and are universal in the sense that they do not depend of the details 
of the model, provided some generic requirements are satisfied (e.g. the variables w have 
only local correlations, and enough non-trivial finite moments, see also the review [23]). The 
models we consider below are predicted to belong to the so-called the Kardar-Parisi-Zhang 
(KPZ) universality class [H] for which ^ = 2/3 and x = 1/3 [MIEB]- In particular, since 
the transition from strong to weak disorder happens at /3 = 0, the universality conjecture 
predicts that for any /3 > the scaling exponents ^, x and the fluctuation statistics equal the 
corresponding zero temperature ones (/3 = oo). 

The first step to this conjecture is to identify x and describe the scaling limit for poly- 
mers and their free energy. This can be done by studying any model believed to be in the 
universality class. The second step is to show universality, i.e., show that the results are 
independent of the chosen model. In this paper we show such universality for any /3 > for 
two particular polymer models - the semi-discrete and continuum directed random polymers. 

Consider now Z^q{t, x) to be the point-to-point partition function of polymers ending at 
X at time r. Denote by Fjg := ■|ln(Z^Q) the free energy and by /^g the free energy limit 
shape (law of large numbers of FJ^q along different velocities). Then, a stronger version of 
the KPZ universality conjecture is to claim that over all models, there exists a unique limit@ 

iTo"' - ^'"nii^^^"'^^)) ■ (1-1) 

One issue is that, since we are in the strong disorder regime, the quenched free energy differs 
from the annealed one (easier to compute) and there is no general way of determining it. The 
conjectural space-time limit is described in [26] where it is called the KPZ renormalization 
fixed point. Information about this fixed point has generally come from studying exactly 

^More precisely, the uniqueness should be up to a spacc-tinic scaling by parameters which can be calcu- 
lated from the microscopic properties of the polymer such as dPo, /3 and the disorder distribution (see, for 
example, [iSllSD] ). 
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solvable models at zero temperature, /3 = oo, such as last passage percolation, TASEP or 
PNG (see the review [23]). For these models, one indeed finds ^ = 2/3, x = 1/3) and 
although the existence of the full space-time limit under these scalings is not known, for fixed 
time r > the limiting spatial process is the Airy2 process [l2l[53]. This extends the results 
of [HIHT] that show that the limiting one-point distribution (fixed (t, x)) is governed by the 
GUE Tracy- Widom distribution [02] • 

Some of the zero temperature models can be reinterpreted as growth models of interfaces 
belonging to the KPZ universality class (where the free energy plays the role of a height 
function). Thus, the universality conjecture says that the height function (properly rescaled) 
converges to the same limit as the polymer free energy. These other models require initial data 
(i.e., height function profile at time zero) and the description of the KPZ fixed point takes 
the limit of this initial data into consideration. The scalings ^ and x do not change, but the 
limiting statistics refiect the initial data (see the reviews [29|[30] and recent experiments [6Tj). 
For random polymers the change in statistics occurs for instance if one looks at point-to-line 
problem instead. We will look at a different way to change the statistics by considering 
point-to-point polymers with boundary perturbations. 

Polymers with boundary perturbations. The role of initial data for growth processes and 
particle systems can be mimicked in the case of polymers by introducing inhomogeneity 
into the quenched disorder so as to encourage the path measure to spend more time on the 
boundary of the underlying path space. For instance consider a polymer with dP^ given by a 
SSRW (either free or point-to-point) in a quenched disorder formed by i.i.d. centered weights 
Wt^xi for X < t] and i.i.d. boundary weights Wt^t with positive mean, for t > 0. As the mean of 
the boundary weights Wt^t increases, the paths which stay along the boundary for a long time 
tend to have a larger Boltzmann weight, and above a critical threshold, the energy-entropy 
competition is dominated by the boundary attraction so that the path spends a macroscopic 
proportion of its time pinned along the boundary. Note that once a path leaves the boundary 
it can not return. This leads to Gaussian scalings and statistics for the free energy and hence 
takes us beyond the basin of attraction for the KPZ renormalization fixed point. 

However, in a scaling window of perturbation strength around the critical value, the bound- 
ary energy is strong enough to modify the free energy statistics and polymer measure, without 
modifying the scaling exponents. Modulating the strength within the window one sees a tran- 
sition from sub to super critical behaviors (in terms of exponents and statistics). Under the 
scaling (11. ip these critically perturbed polymers should converge to the KPZ fixed point with 
a certain type of initial data, which results in different statistics than in the homogeneous 
disorder case. Information about these statistics originated from the analysis of a few solvable 
zero temperature, /3 = oo, models [8l[l0l[l2l[T71[l8l|38]. 

Contributions of this paper. In this work we consider two polymer models with boundary 
perturbations at positive temperature: the O'Connell-Yor semi-discrete directed polymer 
and the continuum directed random polymer (CDRP). Each of these two models are them- 
selves scaling limits of discrete polymers with general weight distributions when the inverse 
temperature f3 scales to simultaneously with the other parameters scaling to infinity. This 
type of scaling is called intermediate disorder scaling as it moves between weak and strong 
disorder (see [2l-ffl[3H[32] ) . This distinguishes these two models as being somewhat universal 
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in their own rights, and hence provides increased motivation to prove the KPZ universahty 
conjecture in their cases. 

The first contribution of this work is to rigorously establish the KPZ universality conjecture 
(for one-point scalings and statistics) for these two positive temperature polymer models with 
boundary perturbations. In particular, we determine the critical perturbation strength and 
show that the phase diagram of one-point scalings and statistics for these polymers match 
that of the solvable zero temperature models. These results significantly expand upon the 
previous proved results for positive temperature polymer models [S1IIS1I2S1EZIISH] • 

Our results rely on an algebraic framework for the exact solvability of these polymers 
which comes from Macdonald processes and their limits and degenerations [16]. One output 
of that work is an exact formula for the Laplace transform of the partition function for 
the semi-discrete polymer with arbitrary boundary perturbation strength. From this, [T6] 
showed KPZ universality for the one-point scalings and statistics for the unperturbed semi- 
discrete polymer at sufficiently low temperatures. The second contribution of this work is to 
develop a variant on the formulas of [TH] for the semi-discrete polymer which are well-adapted 
for rigorous asymptotic analysis at all temperatures (and for boundary perturbations). This 
requires modifications at a fairly high level of the hierarchy of models arising from Macdonald 
processes. 

The third contribution is the discovery and proof of exact formulas for the Laplace trans- 
form of the partition function for the CDRP with boundary perturbations. The semi-discrete 
polymer has an intermediate disorder scaling limit in which it converges to the CDRP [3T|[32]. 
and these formulas are found via rigorous asymptotics of the aforementioned variant on the 
formulas of (TB]. These formulas display a striking similarity (and a clear limit transition) 
to those of the KPZ fixed point with initial data corresponding to the boundary pertur- 
bations [5]. The CDRP free energy solves the Kardar-Parisi-Zhang (KPZ) stochastic PDE 
and the boundary perturbation corresponds to particular class of half Brownian-like initial 
data. The discovered formulas also provide a new class of statistics for this equation, which 
is believed to model the fiuctuations of randomly growing interfaces. 

Let us now explain in detail our results. 

1.2. O'Connell-Yor semi-discrete polymer free energy. The first result is about a 
semi-discrete directed polymer model introduced by O'Connell and Yor |52] . 

Definition 1.1. An up/right path inM.xZ is an increasing path which either proceeds to the 
right or jumps up by one unit. For each sequence < Si < • ■ ■ < stv-i < t we can associate 
an up/right path from (0, 1) to (r, A^) which jumps between the points (sj, i) and (sj, i + 1), 
for i = 1, . . . , N — 1, and is continuous otherwise. Fix a real vector a = (oi, . . . , oat) and let 
B{s) = {Bi{s), . . . ,Bn{s)) for s > be independent standard Brownian motions such that 
Bi has drift at. Let P and E denote the probability measure and expectation with respect to 
the disorder B{-). 

Define the energy of a path cf) to be 



E{4>) = 5i(si) + (fi2(s2) - fi2(si)) + ■ ■ ■ + {BM{r) - Bn{sm-i)) 
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Then the O'Connell-Yor semi- discrete directed polymer partition function Z^(r) is given by 

Z^(r) = J e^^'^UcI), (1.2) 

where the integral is with respect to Lehesgue measure on the Euclidean set of all up/right 
paths (i.e., the simplex of jumping times < Si < ■ ■ • < stv-i < t). The quenched free 
energy F^(r) is defined as 

F^(r) =lnZ^(r), (1.3) 
whereas the annealed free energy is given by F^j(r) := InEZ^(r). 

This model is related to the discrete polymer discussed in the introduction as follows: By 
rotating an angle of 7r/4 the discrete polymer becomes a measure on up/right lattice paths 
starting at (1,1) with weights in (Z>o)^- Consider the point-to-point partition function for 
paths going from (1,1) to (re~^,iV) and rescale the lattice weights like e^^'^Wij. Then as 
e — )■ 0, due to the invariance principle the up/right lattice paths become like (j) of the above 
definition, and the discrete path integral energy becomes E{(j)). 

We focus on scaling the semi-discrete polymer as r = kN for some k > 0. Due to Brownian 
scaling, the parameter k can be changed to 1 by replacing -^(0) by K^^'^E{(j)) in the exponential 
of (ll.2p . In this way, k corresponds to an inverse temperature parameter (3 = k"!'^ . 

It is easy to calculate the annealed free energy exactly as 

Ffn(r) = ^ + (iV-l)lnr-ln(iV-l)! 

which after setting r = kN implies that the annealed free energy density converges to 

f,n(K) := hm iV-iFf„(«:iV) = ^ + ln«: + l. 

Let us briefly relate what the law of large numbers is for the quenched free energy of the 
unperturbed polymer. 

Definition 1.2. The Digamma function is given by ^(z) = ^lnr(2;). For 9 E M+, define 

Ke := vl>'(^), /, := d^>'{d) - ^{6), c, := (-v[/"(^)/2) ^3, 
or, equivalently, for n G M+, define 

r := e M+, r ■■= inf(«:t - ^{t)) = fe^, c'^ := Cg.. 

The following law of large numbers for {nN) was conjectured in [52] and proved in [18]. 
Consider the semi-discrete directed polymer with drift vector a = (0, ... ,0). Then for all 
ft > 0, 

lim N-^¥^{kN) = a.s. 

It follows that for all k > 0, fan('^) = f + lnK+l> f^, and the quenched and annealed 
free energy density converge in the k — ?■ limit. This is in agreement with strong disorder. 

Below we analyze the large asymptotics of the fluctuations of F^{nN) when centered 
by Nf\ 
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Theorem 1.3. Consider the semi- discrete directed polymer with drift vector 
a = (ai, . . . , Qm, 0, . . . , 0) where m < N is fixed and the m non-zero real numbers ai, 02, • . • , ctm 
may depend on N. We can consider without loss of generality that ai> a2> ■ ■ ■ > dm the 
free energy is invariant under permuting these parameters. Then for all k, > 0, we have the 
following characterization of the limiting behavior of the free energy F^{kN) as N ^ 00. 



where -Fbbp,6 is the Baik-Ben Arous-Peche distribution from spiked random matrix 
theory, with b = (61, ... , bm) ■ 

The definitions of -FgueI'") and -Fbbp,6 are provided below in Definition 11.61 The fact that 
the resuh is independent of the ordering of ai, . . . , oat is apparent from the formulas, see e.g. 
Theorem 11.171 below. In Section |2] we reduce the proof of this result to a claim on certain 
asymptotics of the Fredholm determinant formula presented in Section 11.41 We provide a 
formal critical point derivation of these asymptotics in Section [2] and a rigorous proof later 
in Section |5l 

Remark 1.4. If limAr^oo N^^^(ai{N) —9'^) = 00 then the boundary perturbation overwhelms 
the free energy fiuctuations and the scalings and statistics become Gaussian in nature. This 
super critical regime is proved for zero temperature polymers [8l|55]. We do not include a 
proof of this regime here, but it should be readily accessible from the exact formulas via 
asymptotic analysis. 

Remark 1.5. In the unperturbed case, i.e., a= (0,...,0), a tight upper bound on the 
exponent for the free energy fiuctuation scalings was determined in [JS]. In the full one- 
point scaling limit was proved for n > k* with n* a large (enough) constant forced by some 
technical consideration in the asymptotic analysis. In this article we do away with those 
technical issues which allows us to rigorously extend the asymptotics to all positive k,, as well 
as to K tending to zero simultaneously with A^, as we soon will consider. 

Theorem 11.31 is expected by universality, because the same results hold in the zero- 
temperature limit {(3 = 00) and the phase transition is expected to be at /3 = 0. More 
precisely, the limit of the free energy (divided by /3) as /3 goes to infinity and (A^, r) is fixed, 
is described by the ground-state maximization problem 



(a) //limsup^^^iVi/3(ai(iV)-^'') 



—00, then 




where -Fque is the GUE Tracy- Widom distribution I6^ . 
(b) // lim^_,oo N^^^{ai{N) - 9^) = 6^ G M U {-00} for i = 1, . . . ,m, then 





max 

0<si<---<S]v-i<r 



5i(si) + (52(S2) - B2{si)) + ■■■ + {Bm{t) - Bn{sn-i)) . 
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In the unperturbed case, = 0, M^(l) is distributed as the largest eigenvalue of an x 
GUE random matrix [11], see also Theorem 1.1 of [51]. In fact, as a process of r, M^(r) 
has the law of the largest eigenvalue of the standard Dyson Brownian motion on Hermit- 
ian matrices (the lower eigenvalues are also connected to certain generalizations of the free 
energy [51]). It follows from the original work of Tracy and Widom [62] and also [33l[50] that 

For general drift parameters, M'^{t) is related to Dyson Brownian motion with drifts and 
the distribution of {t) then coincides with the largest eigenvalue of a spiked GUE matrix, 
for which the analog of Theorem 11.31 was proved by Peche [55] . The first such results were 
for spiked LUE matrices in the work of Baik-Ben Arous-Peche [5]. 

We now record the definitions of these limiting distributions in terms of Fredholm deter- 
minants. Note that there are many equivalent ways to rewrite these formulas (cf. [7]) and we 
use the most convenient for our purposes. 

Definition 1.6. The GUE Tracy- Widom distribution is defined as 

FcVEir) = det(l - KAi)L2{r,oo), 

where K^i is the Airy kernel, that has integral representations 

KAi{v,v') = JTT^ dw dz ^3/3_w = / c?AAi(r/ + A)Ai(V + A), 

where in the first representation the contours z and w must not intersect. 
The BBP distribution from spiked random matrix theory is defined as 

FBB^,b{r) = det(l - -ft'BBP,6)L2{r,oo), 

where 6 = (&i > &2 > • • • > bm) G and K-BBP,b is given by 

where the path for w passes on the right of bi, . . . ,bm and does not intersect with the path 
z, see Figure [21 It is convenient to extend this definition to allow for bi = — oo for all 
i = i + 1, . . . ,m. Calling b = [bi, . . . ,bi) the finite values of b, we then set i^BBP,6 = -^bbp 1 
and likewise define FBBP,b{r) = F^BPbi''^)- Notice that if m = then A'bbp,;; = -^Ai- For 
representations of this kernel in terms of Airy functions see JEj- 

1.3. Continuum directed random polymer free energy. Just as Brownian motion 
serves as a paradigm for and universal scaling limit of random walks, the continuum di- 
rected random polymer (CDRP) serves a similar role for 1 + 1 dimensional directed random 
polymers [2l[T9]. It is proved that the CDRP is the limit of (general weight distribution) 
discrete [1] or semi-discrete polymers [31] under intermediate disorder scaling in which the 
inverse temperature (or noise strength) is scaled to zero as the system size grows (so as to 
converge to space-time white noise). One such result is quoted as Theorem 13. II 
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Figure 1. Integration contours for w (dashed) and z (solid) of the kernel 
KBBP,b- The black dots are bi, . . . ,b„ 



The directed polymers considered earlier are measures on random walk trajectories 
reweighted by a Boltzmann weight with Hamiltonian given in terms of a path integral through 
a space-time independent noise field along the walk's trajectory. In the continuum, random 
walks are replaced by Brownian motions and the space-time noise field becomes space-time 
Gaussian white-noise. It is possible to define a path measure in the continuum [3]. Here 
we will only focus on the continuum limit of the polymer partition function (11. 2p . not the 
path measure. We define the CDRP partition function with respect to general boundary 
perturbations. These perturbations are the limit of critically tuned boundary perturbations 
of the discrete polymers. 

Definition 1.7. The partition function for the continuum directed random polymer with 
boundary perturbation In^o(^) is given by the solution to the stochastic heat equation with 
multiplicative Gaussian space-time white noise and Zq{X) initial data: 

OtZ = ^dlZ + ZW, Z{0,X) = Zo{X). (1.5) 

The initial data Zq{X) may be random but is assumed to be independent of the Gaussian 
space-time white noise W and is assumed to be almost surely a sigma-finite positive measure. 
Observe that even if Zo{X) is zero in some regions, the stochastic PDF makes sense and 
hence the partition function is well-defined. 

The stochastic heat equation (II. 5p is really short-hand for its integrated (mild) form 

poo pT pco 

Z{T,X)= p{T,X -Y)Zo{Y)dY + / p{T - S, X - Y)Z{S,Y)W{dS, dV) 



where (formally) Gaussian space-time white noise has covariance 



E 



W{T,X)W{S,Y) =6{T-S)6{Y-X), 



and p{T,X) = (27rT)~^/^ exp(— X^/2T) is the standard heat kernel. A detailed description 
of this stochastic PDE and the class of initial data for which it is well-posed can be found 
in [Hlll^, or the review [5^ . 
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As long as is an almost surely sigma-finite positive measure, it follows from work of 
Miiller [H] that, almost surely, Z{T,X) is positive for all T > and X G M. Hence we can 
take its logarithm. 

Definition 1.8. For Zq an almost surely sigma-finite positive measure define the free energy 
for the continuum directed random polymer with boundary perturbation IhZqIX) as 

T{T,X) = \nZ{T,X). 

The random space-time function J-" is also the Hopf-Cole solution to the Kardar-Parisi- 
Zhang equation with initial data J-o(X) = In^o(^) 

The solution to the stochastic heat equation is interpreted as a polymer partition function 
due to a version of the Feynman-Kac representation (see [13] or the review [23]): 

r r /"^ 

Z(T,X) =Eb(t)=x 



Zo(S(0)) : exp : |^ #(t,B(t))di 



where the expectation E is taken over the law of a Brownian motion B which is run backwards 
from time T position X. The : exp : is the Wick exponential which accounts for the fact 
that one can not naively integrate white noise along a Brownian motion (various equivalent 
ways exist to define this exponential [5], [131 [23]). By time reversal we may consider this as 
the partition function for Brownian bridges which can depart at time from any location 
-B(O) e M with an energetic cost of In {p(T,X — B{0))Zo{B{0))) and then must end at X at 
time T. This departure energy cost is the limit of the boundary perturbations for the discrete 
and semi-discrete polymers. 

Let us introduce the class of boundary perturbations which arise in the limit of the semi- 
discrete directed polymer partition function with the first few drift parameters 
tuned as a function of X in a critical way (and all other drifts zero). 

Definition 1.9. The continuum directed random polymer partition function with m-spiked 
boundary perturbation corresponds to choosing initial data for lil.5\) as follows: Fix m > 1 
and a real vector b = {hi, ... , hm); then 

Zo(X) = Z'"(X)lx>o 

where Z'^{X) is defined as in U.S\) with drift vector h, and where lx>o the indicator 
function for X > 0. When m = we will define 0-spiked initial data as corresponding to 
Zq{X) = lx=o, where lx=o is the indicator function that X = 0. 

Theorem 1.10. Fix m > and a real vector b = (6i, . . . , bm)- Consider the free energy of 
the continuum directed random polymer with m-spiked boundary perturbation with drift vector 
b (Definition \1.9\) . 

(a) If m = 0, then for any T > and any S with positive real part, 

(b) If m > 1, then for any T > and any S with positive real part, 
The kernel in the above theorem are given now. 
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Definition 1.11. Fix m > and a real vector b = {bi, . . . ,bm) . The integral kernel 



CDRP,6 



(27ri) 



dw I dz 



an 



^z^/z-zn' ^ Y{aw - bk) 



sin(7r(2; - w)a) e^^l^-'"'^ T{az - bk 



;i.6) 



where a = When m = the product of Gamma function ratios is replaced by 1 and 

the resulting kernel is denoted -RTcdrp- The w contour is from — ^ — ioo to —-^ + ioo and 
crosses the real axis on the right ofbi/a, . . . , bm/cr- The z contour Cz is taken as shifted to 
the right by (see Figure\^for an illustration). Just as in Definition there exist integral 
representations for these kernels involving Airy functions. In particular, 

S 



CDRP 



(r/,r/') 



dt- 



Ai{t + 7])Ai{t + r]'). 



■ S + e-'/- 

Similar formulas exist for KcDKP,b involving Gamma deformed Airy functions f25 \ l3 ^ . 



;i.7) 



c 



4(T 4ct 



Figure 2. Integration contours for w (dashed) and z (solid) of tlie kernel 
KcDKP,b- The black dots are bi/a, . . . , bm/cr. The path equals ^ + C^, so 
as to avoid the zeros of the sine function. 

Remark 1.12. To recover the case of (m — l)-spiked from the m-spiked boundary pertur- 
bation case, one needs to take b^, —oo and simultaneously replace S by —b^S. 

In Section |3l Theorem 13.11 we explain how the semi-discrete polymer partition function 
limits to that of the CDRP. This reduces the proof of the above theorem to a claim on the 
asymptotics of the Fredholm determinant formula presented in Section 11.41 We provide a 
formal critical point derivation of these asymptotics in Section [3] and a rigorous proof later 
in Section |6l 

Remark 1.13. The above Laplace transform can be inverted via a contour integral in S so 
as to give the probability distribution for the free energy. Note that the branch cut in 5'*^^^'^'''^ 
in the integrand in f ll.6p should be taken as the negative real axis. 

For m = the free energy probability distribution was discovered simultaneously and 
independently in both [5l[56] and rigorously proved in [5] via Tracy and Widom's ASEP 
formulas [63]465| [67]: the above Laplace transform formula was soon after (non- rigorously) 
derived from the replica trick approach in [191 [27]. For fn = 1 the free energy probability 
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distribution was discovered and rigorously proved in [25] also via ASEP [66]; the above 
Laplace transform formula was later (non-rigorously) derived from the replica trick approach 
in [39] . The general m > 2 result above is, to our knowledge, new and it is not clear how one 
would derive or prove it from ASEP. 

There are other ways to write the kernel as well as the Fredholm determinant in the 
theorem, as can be seen in the above mentioned citations. 

Remark 1.14. It is not necessary to focus just on the free energy at (T, 0). When m = 0, for 
T fixed, J^(T, X) — lnp(T, X) is a stationary process in X [5] due to the fact that space-time 
white noise is statistically invariant under affine shifts. For m > 1, J-'(T, X) is no longer sta- 
tionary, however, a calculation given in Section 1X2] shows that if we let Zq[X) = Z'"(X)lx>o 
for a shifted drift vector h = {hi + X/T, . . . , 6^ + X/T), then 

Z(r,X) = e-^Z(T,0), 

where Z{T,X) solves the stochastic heat equation with initial data Zq{X). 

A corollary of the above theorem is the large T asymptotics of the free energy fluctuations 
for the CDRP with m-spiked boundary perturbation. 

Corollary 1.15. Fix m > and a real vector b = {bi, . . . , bm)- Consider the free energy of 
the continuum directed random polymer with m-spiked boundary perturbation with drift vector 
ab, with a = (Defimtion\r^. 

(a) If m = 0, then for any r G M, 

T^'J [ (T/2)V3 ^ "J - ^GUE(r), 

where Fque is the GUE Tracy- Widom distribution fU^ (see Definition \1.6\) . 

(b) If m> 1, then for any r G M, 

/ J-(T,0)+T/4! < \ _ ^ . . 

where -Fbbp;6 is the Baik-Ben Arous-Peche IS^ distribution from spiked random matrix 
theory (see Definition M.b]) . 



For m = 0, 1 the above corollary is proved in [5] and [25] (respectively). Given the new 
m > 2 formulas, it is straightforward to prove the full corollary as is done in Section 13.31 

1.4. Fredholm determinant formula for semi-discrete polymer free energy. Theo- 
rems 11.31 and 11.101 are proved via asymptotic analysis of a Fredholm determinant formula for 
the Laplace transform of the O'Connell-Yor semi-discrete directed polymer partition function 
which we now give as Theorem 11.171 The formula is written in terms of a Fredholm deter- 
minant. One of the surprising aspects of the known exactly solvable positive temperature 
directed random polymers is that the Laplace transform of their partition functions are given 
by Fredholm determinants. Here we write the Laplace transform as a double-exponential 
transform of the free energy. 
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R 



Figure 3. Left: The contour Cq,^^ (dashed), where the black dots are 
ai, . . . , Oat and a > max{ai, . . . , aj\f}. The contour v + is the sohd hne. 
Right: The contour Vy, where the thick part is T>y^^ and the thin part is T>y^\. 
The grey dots are at {1, 2, . . .} 



Definition 1.16. For a G M and (f G (0,7r/4), we define a contour C a, that surrounds the 
portion of the real axis with values less than a hyCa^^ = {a+e''''^'^'^-'?/}ygK^U{a+e'*^'^^'^''?/}ygiR^. 
The contour is oriented so as to have increasing imaginary part. 

For every v G C^,^ we choose R = — Re{v) + a + 1, d > 0, and define a contour as 
follows: Vi, goes by straight lines from R — ioo, to R — id, to 1/2 — id, to 1/2 + id, to R + id, 
to R + ioo. The parameter d is taken small enough so that v + do not intersect Ca^ip. We 
also call the portion of with real part R and V^^^ the remaining part. See Figure\^ 
for an illustration. 



Theorem 1.17. For fixed N > 1, r > 0, a drift vector a = (ai, . . . , Cat), and 
a > max{ai, . . . , Ctv}, the Laplace transform of the partition function for the O'Connell-Yor 
semi-discrete directed polymer with drift a is given by the Fredholm determinant formula 



E 



:tZ^(T) 



det(l + i^„)i2 



where Ca,ip is given in Definition \1.16\ for any ip G (0,7r/4). The operator Ku is defined in 
terms of its integral kernel 

K^^y^y') = ^ dsTi-s)T{l + s) n -. (1.8) 

Theorem 5.2.10 of [12] gives a similar formula for the Laplace transform of the semi-discrete 
directed polymer partition function. The difference between the two formulas is the contours 
involved - both the contour of the space and the contour involved in defining the kernel. 
The contours of the above formula are unbounded. This somewhat technical modification 
is important since it enables us to perform rigorous steepest descent analysis as necessary 
to prove Theorems 11.31 and 11.101 In [16] , the corresponding contours were bounded, thus 
hmiting the asymptotic analysis of the semi-discrete polymer to a low temperature regime. 
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In order to prove this modified formula we modify the choice of contours very early in 
the proof of Theorem 5.2.10 of [16]. In Section |4] we provide a proof of the above formula 
(see in particular the end of Section 14. ip with the more detailed and technical pieces of the 
proof delayed until Section [71 The occurrence of unbounded contours introduce some new 
considerations in proving this theorem. It is not presently clear how to derive the above 
theorem directly from Theorem 5.2.10 of |16] . 
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Research through the Schramm Memorial Fellowship. PF was supported by the German 
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2. Free energy fluctuations for the O'Connell-Yor semi-discrete polymer 



In this section we reduce the proof of Theorems 11.31 to a statement about the asymptotics 
of a Fredholm determinant (Theorem 12.11 below). We then provide a formal critical point 
derivation of the asymptotics, delaying the rigorous proof until Section O 

The starting point for the proof of Theorem ll.3l is the Fredholm determinant formula given 
in Theorem 1 1 . 1 71 for E[e~"'^ ^*'*]. We rely on the fact that under the scalings we consider, the 
Laplace transform of the partition function converges to the asymptotic probability distribu- 
tion of the free energy (a similar approach is used in the proof of Corollary I1.15p . Towards 
this aim, define a sequence of functions {67v}Ar>i by 6Ar(x) = exp (— exp (c'^A^/^x)), where 

is given in Definition 11.21 Recall also that we are scaling r = kA for k > fixed. 

Assume that the drift vector a = (ai, . . . ,a7v) is as specified in the statement of Theo- 
rem [T3i Set 

where is as in Definition 11.21 and observe that 



E 



E 



^ , F^(kA) - AP 



(2,1) 



The A — oo asymptotics of the left-hand side of (12. ip can be computed from taking asymp- 
totics of the Fredholm determinant formula given in Theorem 11.171 which states that 



E 



-uZ^(KAf) 



det(l + A„)i2 



(2.2) 



for a > max{ai, . . . , qn} and G (0, 7r/4). This asymptotic result is stated as Theorem 12.11 
below and proved in Section [51 After explaining how it implies Theorem 11.31 we provide a 
formal critical point derivation of the asymptotics. 
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Theorem 2.1. Consider a vector a = (ai, . . . , a^, 0, . . . , 0) where m < N is fixed and the 
m non-zero real numbers ai, 02, • • • , "^o?/ depend on N. We can consider without loss of 
generality that Oi > 02 > . . . > a„i > 0. Then for all n > 0, 

(a) For the unperturbed case, m = 0, 

lim det(l + Ku)l^Cc ^) = det(l - /'^Ai)L2(r,oo) = ^GUE(r-), 

where -Fgue is the GUE Tracy- Widom distribution I6^ . 

(b) // lim^^oo N^/'^^iaiiN) - ^'^) = 6^ e M U {-00} for i = 1, . . . ,m, then 

lim det(l + Ku)l^Cc, ^) = det(l - i^BBP,b)L2(r,oo) = FBBP,b{r), 

where Fbbp,6 is the Baik-Ben Arous-Peche |^ distribution from spiked random matrix 
theory. 

We remind from Definition 11.61 that when hi = —00, 1 < i < m, then -Fbbp,6 = -^gue- 
The above resuh imphes the the right-hand side of (12. ip has a hmit p{r) that is a continuous 
probabihty distribution function. Here p{r) is the hmiting distribution function in cases (a) 
or (b) of Theorem 12.11 The functions QnIx — t) approximate l{x < r) in the sense necessary 
to apply Lemma 18. Ij and hence p{r) also describes the limiting probability distribution of 
the free energy: 

'F^(KAr) - Ar/« 

nm IT I 

This implies Theorem 11.31 



lim P -— < r = p(r) 



2.1. Formal critical point asymptotics for Theorem 12.11 We provide a formal analysis 
of the asymptotics of the Fredholm determinant det(l + Ku)L'^{(Za^)- particular, we only 
focus on the limit of the kernel Ku, and even in that pursuit, we only consider the pointwise 
limit of the kernel in (11. 8p . We also disregard issues respecting the choice of contours. All of 
these issues are considered in the rigorous proof contained in Section [5l 

The first set of manipulations to that we make are to rewrite r(— s)r(l + s) = 
— tt/ sin(7rs), and to factor out the ratios of Gamma functions involving those for 1 < i < m. 
With a change of variable 5 = s + t', we obtain 



-1 f vr exp (iVG(.) + rc-iVV3,) ^ ^ r(. - a.)r(.-) 



2m J sm{7i{z - v)) exp {NG{z) + rc'^N^/^z) z - v' T{z - ak)T{v] 
where 



(2.3) 



G{z)=\nT{z)-Kj + f^z. (2.4) 

The problem is now prime for steepest descent analysis of the integral defining the kernel 
above. The idea of steepest descent is to find critical points for the function in the exponential, 
and then to deform contours so as to go close to the critical point. The contours should 
be engineered so that away from the critical point, the real part of the function G in the 
exponential decays and hence as N gets large, has negligible contribution. This then justifies 
localizing and rescaling the integration around the critical point. The order of the first non- 
zero derivative (here third order) determines the rescaling in N (here N^^^) which in turn 
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corresponds with the scale of the fluctuations in the problem we are solving. It is exactly this 
third order nature that accounts for the emergence of Airy functions and hence the Tracy 
Widom (GUE) distribution as well as the BBP transition distributions. 
The critical point equation for G is given by G\z) = with 

G'{z) = ^(z) - Kz + r. 

The Digamma function "^{z) = ^ln(r(z)) is given in Definition 11.21 Also given in that 
definition is 6'^ G M.^ which is the critical point, i.e., G'{6'^) = 0. At the critical point 
G'^e"") = and G^^^e"") = ^"{9'') = -2{c'^f so that Taylor expansion at the critical point 
gives (up to higher order terms) 

G{v) - G{9^) - ^(i; - e^Y, G{z) ~ G{e^) - - 9^f. 

3 3 

This cubic behavior suggests rescaling around 6'^ by the change of variables 

Under the above change of variables we find that as — )■ oo, 

exp {NG{v) + rc'^N^'^v) exp (z^S - rz) 
e'x^{NG{z) + rc'^N^/^z) exp (ty'^/S — rii?) 

Note that since the v,v' variables were scaled, there is a Jacobian factor of l/(c'^A^^/'^) in- 
troduced into the kernel. Grouping this with the reciprocal sine function we see that as 
N ^ oo, 

1 TT 1 dz dz 

c'^A^i/s sin(7r(f — ^)) w — z^ z — v' z — w' 

It remains to study the ratio of Gamma functions. This is where the subcriticality (a) versus 
criticality (b) becomes important. First notice that the factor YYk=i ~^ ^- The fact 

that the critical value for the a^'s is 6'^ coincides with the centering of the change of variables 
is not an accident as we now explain. After the above change of variables 

Tiv - ak) _ T{e^ -ak + w/{c''N^/^)) 

r{z-ak) ~ r{e- - + z/{c'^m/^)) ' 

As long as lim supjv_>.oo N^^^i^-kiN) —6'^) = — oo, as — > oo, the numerator and denominator 
both converge to r{6'^ — ak) and hence their ratio is 1. Thus for the subcritical case (a) the 
limiting kernel is given by 

1 1 exp (2^/3 — r^;) 

27ri J {w — z){z — w') exp (w^/S — rw) 

In the critical case, limsup^_j.oQ A^"'^/^(afc(A^) — 6'^) = bk ioi 1 < k < m. This means that 
after the change of variables 

T{v - ak) _ T{{w - h)/{c^N'/^)) _ z - 
T{z - ak) T{{z - bk)/{c-m/^)) ~ w-bk 
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for large A^, since T{z) ~ l/z near 0. Therefore, in the critical case (b) the limiting kernel is 
given by 

~ , 1 /■ 1 exp{z^/3-rz) z -bk 

KBBP,b{'w,w) = — / dz- — — -M -. 

zTTi J [w — z)[z — w') exp — rw) w — Ok 

The variables w,w' are integrated along a contour C from e~^'^'/^oo to e'^'^^^'^oo, passing on 
the right of 6i, . . . , b^- The variable z is integrated along a contour which goes from e~'^^^^oo 
to e'^'/^oo without crossing C. 

The limiting Fredholm determinants we have derived can be rewritten as Fredholm deter- 
minants on L^(r, oo) as shown in Section I8.3[ and hence one sees their equivalence to those 
given in Definition 11.61 

This completes the formal critical point derivation of Theorem 12.11 Let us note, that in 
order to make this formal manipulations into a proof it is necessary (among other things) to 
be careful about the contours. First one has to find a steep descent pathQ for G{v), which 
might not be obvious due to the Digamma function (it turns out a useful representation 
for the real part of the Digamma function is as an infinite sum, see Section I3TT]) . Secondly, 
one would like to find a steep descent path for —G{z), but because the path for z has to 
include all the poles at f + l,f + 2,..., such path does not exists. The way out we used was 
to find a steep descent path for —G{z) and then add the contributions of the poles at the 
V + 1, . . . ,v + i which lie on the left of the path (see Figures |6] and [7]). Finally, one needs 
to get estimates so that not only the kernel converges, but also the Fredholm determinant. 
Further technical details are presented in the proof, see Section [51 



3. Laplace transform of the CDRP partition function 

In this section we reduce the proof of Theorem II. 101 to a statement about the asymptotics 
of a Fredholm determinant (Theorem 12.11 below). We then provide a formal critical point 
derivation of the asymptotics, delaying the rigorous proof until Section [61 We also include 
two brief calculations delayed from the introduction. 

The CDRP occurs as limits of discrete and semi-discrete polymers under what has been 
called intermediate disorder scaling. This means that the inverse temperature should be 
scaled to zero in a critical way as the system size scales up. For the discrete directed polymer 
it was observe independently by Calabrese, Le Doussal and Rosso [19] and by Alberts, Khanin 
and Quastel [2] that if one scaled the inverse temperature to zero at the right rate while 
diffusively scaling time and space, then the discrete polymer partition function converges 
to CDRP partition function. Using convergence of discrete to continuum chaos series, |1] 
provide a proof of this result that is universal with respect to the underlying i.i.d. random 
variables which form the random media (subject to certain moment conditions). 

Concerning the semi-discrete directed polymer, [21] prove convergence of the partition 
function to that of the CDRP. The results of [21] deal with zero drift vector, and in [32] the 

■^For an integral I = dze*^^''\ we say that 7 is a steep descent path if (1) Re(/(z)) reaches the maximum 
at some zq g 7: Re(/(2:)) < Re(/(zo)) for z <S 7 \ {zq}, and (2) Re(/(2)) is monotone along 7 except at its 
maximum point zq and, if 7 is closed, at a point zi where the minimum of Re(/) is reached. 
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convergence is extended to deal with a finite number of non-zero drifts, critically tuned so as 
to result in an boundary perturbation for the CDRP partition function. 

Theorem 3.1 ( |3lll32]). Fix T > 0, X e R, m > 0, and a real vector b = (fei, . . . 
Set K = ^/T/N from which t = kN = VTN, and 6 = 6'^^ \/N/T + \. For each N >m 
define a drift vector a = (oi, . . . , a^, 0, . . . , 0) where ak := 9 + for 1 < k < m. Consider 
the O'Connell-Yor semi- discrete polymer partition function Z^(r) with drift vector a. For 
yTN + X > 0, define its resettling as 

^ _ z^ivm+x) 

^ C{N,T,X) 

with scaling constant 

C{N,T,X) = exp (^N + lNln{T/N) + ^{VtN + X) + X^nJt^ exp {-\mhi{T/N)) . 

Then as N oc, Z^{T^X) converges in distribution to Z{T,X) which is the partition 
function for the continuum directed random polymer with m-spiked boundary perturbation 
corresponding to drift vector b. 

Remark 3.2. For S with positive real part, due to the almost sure positivity of Z{T,X), 
the above weak convergence also implies convergence of Laplace transforms: 

E re-^^(^'^)l = hm E [e-^^"(^'^)" 

Let us focus on the semi-discrete polymer with drift vector as specified in the above theo- 
rem, and X = 0. Then, rewrite the CDRP partition function in terms of the free energy, the 
previous remark implies 

where 

^ ^ ^g-7V-|7Vln(T/7V)-iv^+T/4!g|mln(T/iV)^ ^^^^^ 

a > max{ai, . . . ,aAr} and ip G (0,7r/4). This reduces the proof of Theorem 11.101 to the 
following. 

Theorem 3.3. Fix S with positive real part, T > 0, m > 0, and a real vector b = [bi, . . . , b^) . 
Set K = ^JT/N from which r = kN = VTN, and 6 = 6''^ ^JN/T + \. For each N > m 
define a drift vector a = (ai, . . . , am, 0, . . . , 0) where ak := + bk for 1 < k < m, and define 
u by ^3. ij) . Fix any a > max{ai, . . . ,aAr} and (f G (0,7r/4). Then, recalling the kernel Ku 
from U.S\) . it holds that: 

(a) In the unperturbed case, m = 0, 

lim det(l + Ku)l^c^ ^) = det(l - Kcdrp)l2{r+) 



E [e-^-P(-^(^'0)+r/4!)] ^ ^lij^E ^g-«z-(VT]v)j ^ Jim det(l + K^)lhc^^^) 



with -K'cDRP given in Definition 
(b) For the perturbed case, m > 1, 



lim det(l + Ku)l^{c^ \ = det(l - /'^cDRP,b)L2(R+) 



with -R'cDRP b given in Definition 
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This theorem is proved in Section |6l However, we now include a formal critical point 
derivation of the asymptotics. 



3.1. Formal critical point aysmptotics for Theorem 13.31 We provide a formal analysis 
of the asymptotics of the Fredholm determinant det(l + Ku)L^{Ca^f) under the prescribed 
scalings. In particular, we only focus on the limit of the kernel Ku, and even in that pursuit, 
we only consider the pointwise limit of the kernel in (11. 8p . We also disregard issues respecting 
the choice of contours. All of these issues are considered in the rigorous proof contained in 
Section [61 

Set K = ^/tJN and note that u = (5/0™)e-^/''+^(^"'^') (recall Definition O]) so that 
the result is the same if we just set u = {S/9^)e~^^'^ . It is convenient to do the change 



of variable v ^ 9 + aw and z — 
approximations) given by 

-1 /• , OTXS^' 

/ dz — — — — 

27ri / sin(7r(2 



9 + oz. In these new variables, the kernel is (up to the 



Koiw, w' 



-w)a NG{e+c7w)-NG{e+c7z) ™ "p 

-n 



aw 



az] 



wja 



w 



k=l 



T{az - bk)T{9 + aw)9" 



where G is given by (12. 4p . 

G has a double critical point at 9, as G'{9) 
Taylor expansion of G one sees that 



G"{9) = 0. Therefore, using the large N 



NG{9 + aw) ~ NG{9) 



-w 



Also, as 9 is going to infinity with A^, it is immediate that r( 
Thus, as A^ — 7- oo, the kernel Kq goes to 



az)9'''^ /T{9 + aw)9''' 1. 



A'cDRP,fe(w^, w') 



-1 
27ri 



dz 



an 



g{z-w)a g.3/3_^3/3 m ^^^^ _ 



sin(7r(z — w)a) z — w' 



n 

k=l 



Tiaz - hk) 



The variables w, w' are on the contour and z on of Figure [2l The Fredholm determinant 
with this kernel can be rewritten as a Fredholm determinant on L^(M+) as shown in Section [873] 
and one checks that this corresponds with det(l — A'cDRP,b) l2(m+) as desired. 

This completes the brief and formal critical point derivation of Theorem [3731 The technical 
challenges in order to produce a rigorous proof are similar to those explained at the end of 
Section [2l There is a new difficulty with regards to contours which arises in this case, 
however. The variables w, w', z arose from the change of variables of w', z. Besides a shift 
by 9 (which goes to infinity when A^), the variables are just scaled by a, which is of order 1 
(here T is fixed). This means that all of the singularities of 1/ sin(7r(5 — v)) remain in the 
A^ — oo scaling. This necessitates significant care in the choice of contours along which to 
take asymptotics. The complete proof of these asymptotics is given in Section [61 

3.2. AfRne shifting of the perturbed CDRP partition function. The following calcu- 
lation shows how the claim (for m > 1) in Remark [1.1 41 is derived. As explained in Section [T73l 
we can write 



Z{T,X) 



p(T,F-X)Zo(F)EB(o)=y 

) B{T)=X 



exp : <^ I W{t,B{t))dt 



dY, 
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where Kb{o)=y denotes the expectation over a Brownian bridge B starting at -B(O) = Y and 

B{T)=X 

ending at B(T) = X. Let B{s) = B{s) — sX/T, then observe 



Z{T,X) 



piT,Y) 



p{T,Y) 



p{T,Y - X)ZoiY)E, 



exp 



W{t,B{t))dt 



dY, 



■'B{0)=Y 
B{T)=0 

where the space-time white noise here is the affine shift of the previous one (but, in any case, 
equal in law). We have also inserted a factor of 1 so that we can now rewrite 

r-T 



Z{T,X) 



e 



p(T,F)Zo(F)e^^/^E 



B(o)=y 

i3(T)=0 



exp 



W{t,B{t))dt 



dY. 



If Zo{Y) = Z™(F)ly>o for drift vector b = (61,..., 6„) then Zo{Y) = Zo{Y)e^^/^ = 
Z™(F)ly>o for drift vector b = {bi + X/T, . . . , 6^ + X/T). Letting Z correspond to the 
solution to the stochastic heat equation with Zq{Y) initial data, we find that 



z(r,x) 



e-^Z{T,0). 



Thus, after a parabolic shift, and an addition of drift into the boundary perturbation, the 
distribution of the general X free energy can also be determined from Theorem ILIOI as well. 

3.3. Proof of Corollary 11.151 Recall a = Let us focus on m > 1 and note that 

m = is proved identically. Define a sequence of function {0t}t>o by 63- (x) = e'^""'" . Then 
if we set S = e~'^l" 



E {t 



-5exp(J-(T,0)+T/4!) 



] = E 



e 



T 



TiT, 0) + T/4! 



(7 



-1 



(3.2) 



Let us first calculating the T — > 00 limit of the Fredholm determinant expression for the 
left-hand side of the above equality. It is easy to see that in this limit 

-fi^CDRP,<7fe(?7, V) -> -K'bbp,6(?7 + r, 77' r). 
This is because T — j- 00 corresponds to a — !■ and thus 

^^Q{z-w)a ^-r{z-w) Y{a{w - 6fe)) ^ z - bk 



sm TTU — w]a] 



z 



w 



T a z 



bk)) ' w-bk 

One readily identifies the resulting expression with that of Definition 11.61 The tail bounds 
necessary to justify this are not hard (see [5l[25] for instance). Going back to (13. 2p . this 
implies that 

[F{T, 0) +T/4! 



lim E 

T-5.00 



FBBY'.bi.n- 



(3.3) 



Since (away from x = r) Qxix — r) is converging to lx<r (and in particular due to Lemma IHTTll . 



lim E 

T-s>oo 



6 



T 



fJ^{T,0) + T/A\ 



a 



lim P 

T-s>oo 



JP(T, 0) + T/4! 



< r 



and hence the corollary is proved. 
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4. Proof of Theorem I1.17[ the main formula 



The proof of Theorem II .171 follows closely the proof of Theorem 5.2.10 in [16]. The major 
difference is that the formula resulting from Theorem 5.2.10 had bounded contours which 
were unsuitable for the full scope of asymptotic analysis necessary to prove Theorems 11.31 
and 11.101 This somewhat minor modification to the final formula requires us to modify the 
proof quite early on and in fact the unboundedness of contours results in a fair number of 
new technical steps in the proof. We produce in this section the complete proof. Results 
used along the way which are stated and proved in [TB] are not reproved. Some of the more 
technical points in the proof are delayed until Section [3 

To prove Theorem 11.171 we use the theory of Macdonald processes as developed in [T6] . 
As we explain below in Section 14. H due to O'Connell's work [51] on a continuum version 
of tropical RSK correspondence, the partition function Z^(t) arises as a marginal of the 
Whittaker process (or measure). Macdonald processes sit above Whittaker processes due to 
the hierarchy of symmetric functions. Under suitable scaling, Macdonald processes converge 
weakly to Whittaker processes (and hence a suitable marginal converges to Z^(t)). Due to 
the Macdonald difference operators and the Macdonald version of the Cauchy identity it is 
possible to compute simple formulas for expectations of a large class of observables of Mac- 
donald processes. In particular it is possible to compute g-moments for the marginal random 
variable which converges to Z^(t). A g-Laplace transform can be rigorously computed by 
taking an appropriate generating function of the g-moments, and switching the expectations 
and summation is rigorously justified. For the g-Laplace transform we find a nice Fredholm 
determinant. Taking the degeneration of Macdonald processes to Whittaker processes, the 
g-Laplace transform becomes the usual Laplace transform, and due to the weak convergence, 
we recover the desired formula for the Laplace transform of Z^(t). 

We do not provide an introduction to the theory of Macdonald symmetric functions here. 
Instead we refer readers to Section 2.1 of [16] for all of the relevant details, or to Chapter 
VI of [17]. We also do not define Macdonald processes in their full generality but content 
ourselves with studying a certain set of marginals of the processes which are called Macdonald 
measures, and a certain Plancherel specialization. 

4.1. O'Connell's Whittaker measure and its relation to polymers. In order to state 
the pre-asymptotic Laplace transform formula and set up its derivation, it is useful to intro- 
duce a few concepts. Initially they may seem a little out of place, but due to O'Connell's 
work [51] (some of which is recorded in Theorem 14. 3p the connection to the semi-discrete 
polymer becomes clear. This connection is analogous to the relation between last passage 
percolation and the Schur process (see e.g. [12] )• 

The class-one qIj^- Whittaker functions are basic objects of representation theory and in- 
tegrable systems [281 HS]. Oiis of their properties is that they are eigenf unctions for the 
quantum g[^-Toda chain. As showed by Givental [36], they can also be defined via the 
following integral representation 




N-l k 
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where A = (Ai, . . . , Xn) and 

N / k 



Tx{x) = Afc ^Xk-i 



k-l 




N-1 k 





k=l \i=l 

For any r > set 



i=l 



k=l i=l 



j 



N 



9riXi, . . .,Xn) 



ijy{xi, XN)e "^^^rriNii^) Y\_ 



with the Skylanin measure 



mAr(l/) 




Note that our definition of Whittaker functions differs by factors of i from those considered 
by O'Connell [S]. 

Definition 4.1. For N > 1, define the Whittaker measure with respect to a vector 
a = {tti, . . . ,a]\f) G via the density function (with respect to Lebesgue) given by 



The fact that this measure integrates to one follows from analytic continuation of the 
orthogonality relation for Whittaker functions (see [16] Proposition 4.1.17). We write expec- 
tations with respect to Whittaker measures as (■)wM(^ ^ 



Let us remark that in [T6l[5T] the Whittaker measure is a level marginal of a measure on 
triangular arrays called the Whittaker process, which is also defined via Whittaker functions, 
but whose precise definition will not be important for us in this work. 

The connection between the Whittaker measure and the semi-discrete polymer is best 
explained by introducing an extension to the polymer partition function and free energy. 



Definition 4.2. The hierarchy of partition functions (r) for < n < N is defined so that 
Zq (r) = 1 and for n > 1, 



where the integral is with respect to the Lebesgue measure on the Euclidean set (r) of all 
n-tuples of non-intersecting (disjoint) up/right paths with initial points (0, 1), . . . , (0,n) and 
endpoints (r, — n + 1) , . . . , (r, A^) . For a path 0j which starts at (0, i) , ends at {t, N — n + i) 
and jumps between levels j and j + 1 at times Sj, the energy E{(j)i) is defined as 



with Bi, . . . ,Bn independent Brownian motions. It follows that Zf^(r) = Z^(r) as defined 
in [L^] . 

The hierarchy of free energies (r) for 1 < n < N is defined via 



(4.1) 
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It follows that Ff (r) = F^(r) as defined m / fO) . 



The following result is shown in [51] by utilizing a continuous version of the tropical RSK 
correspondence (see also [23] for a discrete analog) and certain Markov kernel intertwining 
relations. 

Theorem 4.3. Fix N > 1, r > and a vector of drifts a = («!,..., aj^), then {F^ ('^)}i<n<Ar 
is distributed according to the Whittaker measure WM(_(i^^...^_(j^.t-) of ( [/^. 

More is shown in [51] including the fact that the collection of free energies evolves as a 
Markov diffusion with infinitesimal generator given in terms of the quantum gij^ Toda lattice 
Hamiltonian; as well as the fact that the entire triangular array F(r) is distributed according 
to the Whittaker process which we briefly mentioned earlier. 

Remark 4.4. It is useful to note the following symmetry: The transformation Tj -H- —Tjsf-^-i-i 
maps WWl(^_ai,...-aN;T) to WM(ai ,...,ajv;r) (the sigu of tt/s chauges). This easily follows from 
the definition of WM(ai,,,,,ajv;r)- 

The result below is similar to the formula found in Theorem 4.1.40 of [16] except that the 
contours have been modified. This technical change requires us to go through the proof of 
this theorem and perform a number of estimates which are harder than in [IS], due to the 
unboundedness of the contours. 

Theorem 4.5. Fix N > 1, r > and a vector a = (ai, . . . , a^) G and a > max{aj}. 

Then for all u E C with positive real part 

where the operator Ku is defined in terms of its integral kernel 

K^^y^y') = ^ dsT{-s)Til + s) n -. (4.2) 



The paths Ca,ip and T>^ are as in Definition \l.l(A where ip is any angle in (0,7r/4). 

This theorem is proved over the course of the rest of this section, with the more technical 
aspects of the proof relegated to the later Section [7l Before going into this, let us note that 
in view Theorem 14. 3^ the proof of Theorem 11.171 is now immediate. 

Proof of Theorem \1.17\ By appealing to Theorem 14.31 and Remark 14.41 we relate Z^(r) to T/v 

as 

-«Z^(r) 



e 



Applying Theorem 14.51 gives the claimed Laplace transform and completes the proof. □ 
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4.2. Macdonald measures. We now turn to the proof of Theorem 14. 5[ Before giving the 
proof (as we do in Section 14. 5p we need to recall and develop a variation on the solvability 
framework of Macdonald measures [16] . 

The Macdonald measure is defined as (following the notation of [16] where it is introduced 
and studied) 

MM(ai, . . . , aTv; p)(A) = — -— z ^ • 

n(ai, ...,aN;p> 

Here A is a partition of length < N (i.e., A = (Ai > A2 > ■ ■ • Aat > 0)) and Pa and Qx are 
Macdonald symmetric functions which are defined with respect to two additional parameter^ 
q,t & [0, 1). The notation Px{ai, . . . ,0,^) means to specialize the symmetric function to an 
N variable symmetric polynomial, and then evaluate those variables at the di The 
notation Q\{p) means to apply a specialization p to Qx. A specialization of a symmetric 
function / in the space Sym of symmetric functions of an infinite number of indeterminants 
is an algebra homomorphism p : Sym — )■ C, and its application to / is written as /(p). We 
will focus here on a single class of Plancherel specializations p which are defined with respect 
to a parameter 7 > via the relation 

J^dnipK = exp(7M) =: Il{u;p). (4.3) 

n>0 

Here u is a formal variable and Qn = Q{n) is the (g, t)-analog of the complete homogeneous 
symmetric function hn- The above generating function in u is denoted H(u; p). Since Qn forms 
a Q[q',t] algebraic basis of Sym, this uniquely defines the specialization p. The Plancherel 
specialization has the property that Qx{p) > for all partitions A. 
On account of the Cauchy identity for Macdonald polynomials 

^ Px{di, . . . , dN)Qxip) =■■ n(ai, . . . , Stv; p) = H(ai; p) ■ ■ ■ n(a7v; p). 

\:e{X)<N 

Our ability to compute expectations of observables for Macdonald measures comes from 
the following idea: Assume we have a linear operator T> on the space of functions in 
variables whose restriction to the space of symmetric polynomials diagonalizes in the basis 
of Macdonald polynomials: VPx = dxPx for any partition A with £(A) < A^. Then we can 
apply T) to both sides of the identity 

P\{ai, dN)Qx{p) = n(ai, . . . , Oat; p). 

A:^(A)<Ar 

Dividing the result by H(ai, . . . , oat; p) we obtain 

, , \ PH(5i,...,ajv;p) 

\aA;MM(ai,...,ajv;p) " -^0^ ~ T ' '^^■'^^ 

i.i.[ai, . . . , a^; p) 



^The parameter t does not correspond to time in any of the processes we have considered. For time we 
have either used r in the semi-discrete polymer or T in the CDRP 

^We use tildes presently since after a limit transition of Macdonald measures to Whittaker measures the 
Si's will become a^'s 
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where (■)MM(ai,...,a]v;p) represents averaging ■ over the specified Macdonald measure. If we 
apply T) several times we obtain 

. _ P^n(ai,...,ajv;p) 

Ill^ai, . . . , Cat, p) 

If we have several possibilities for T) we can obtain formulas for averages of the observables 
equal to products of powers of the corresponding eigenvalues. One of the remarkable features 
of Macdonald polynomials is that there exists a large family of such operators for which they 
form the eigenbasis (and this fact can be used to define the polynomials). These are the 
Macdonald difference operators. We need only consider the first of these operators D\i which 
acts on functions f oi N independent variables as 

. , .,. Xi Xj 

This difference operator takes symmetric polynomials to symmetric polynomials and acts 
diagonally on Macdonald polynomials. 

Proposition 4.6 (VI(4.15) of [17|). For any 'partition A = (Ai > A2 > ■ ■ ■ ) with A^ = /or 

D\Px{x^, ■■■ ,x^) = + g"^t^-2 + ■ ■ ■ + g^-) P,(a:i, ■ ■ ■ , x^). 

Although the operator D\q does not look particularly simple, it can be represented by 
contour integrals by properly encoding the shifts in terms of residues. 

Proposition 4.7 (Proposition 2.2.13 of [H]). Fix k > 1. Assume that F[ui, . . . ,un) = 
f{ui) ■ ■ ■ f{ui\f). Take Xi, . . . ,X]sf > and assume that f{u) holomorphic and nonzero in a 
complex neighborhood of an interval in M that contains {q'^Xj 
Then 

({DifF){x) _ (t-l)-" r r (tz^-qzt){z^-zt) -X, ^1 

Fix) (2.i)- ; ■ ■ JJl^^ (z. - cz,) (tz.- z,) 11 1^11 

where the Zc-contour contains {qzc+ii ■ ■ ■ iQZkiXi, . . . ,xn} o,nd no other singularities for 
c = 1, . . . ,k. 

Observe that n(ai, . . . , qn] p) = Hill n(aj; p). Hence, Proposition 14.71 is suitable for eval- 
uating the right-hand side of equation (14. 4p and hence computing the associated observable 
of the Macdonald process. 

4.3. The emergence of a Fredholm determinant. Macdonald polynomials in variables 
with t = are also known of as q-deformed gljy Whittaker functions [35j. We now denote the 
Macdonald measure as MMf=o(Si, . . . , a^', p) and refer to these as q-Whittaker measures. 

The partition function for the corresponding g- Whittaker measure MMf=o(ai5 ■ ■ ■ ,5,n] p) 
when p is Plancherel simplifies as 

N 

Px{ai, . . .,aN)Qx{p) = n(ai, ...,dN]p) = ]^exp(7aj). 

AeY(Af) j=l 
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4.3.1. Formulas for q-moments. Let us take the limit t — ?• of Proposition 14.71 Write 



Then 



(o^,^k f - f 11 ....... dzi---dzk, (4.5) 



(27ri)''' J J f-^ Za — qzs zi ■ ■ ■ Zk 



where f{z) = YliLi 5^ 6xp{(g— 1)7,2} and where the Zj-contours contain {qzj^i, . . . , qzk} as 
well as {Si, . . . , a^r} but not 0. For example when k = 2 and all di = 1, Z2 can be integrated 
along a small contours around 1, and Zi is integrated around a slightly larger contour which 
includes 1 and the image of q times the Z2 contour. 

We have encountered the point at which we diverge from [16]. In particular, we will 
now deform our contours from bounded curves to infinite contours (as justified by Cauchy's 
theorem). This may not seem so significant presently, but in the final formula this frees 
us up to deform contours to steepest descent contours and hence to rigorously prove the 
various asymptotics desired for our limit theorems. However, the unboundedness of contours 
introduces new complications which we must address. 

First, however, let us define two important unbounded contours we will soon encounter. 

Definition 4.8. For an illustration, see Figure^ For any a > and (f G (0,7r/4) define 
C^^^ = {a + e-''^'f5n(s/)^^^ g ]^} ^^oie that it is oriented so as to have decreasing imaginary 
part). For w G Ca,ip for some choice of parameters, define a contour in the same way as 
Vyj of Definition \1.16\ hut with R and d taken such that the following holds: For all s G Vyj 
(i) if b = J — |, then aTg{w{q^ — 1)) G (7r/2 + 6, 37r/2 — b); (ii) q^w lies to the left ofCa,^- 

Remark 4.9. Let us check that the contours V^, in above definition actually exist. Fix 
f G (0,7r/4) and fix a contour Ca,^, as above. For w G Cq,_<^ we must show that there exists 
R and d as desired to satisfy (i) and (ii). The existence of these contains is clear for \w\ 
small, so let us consider when |t<;| is sufficiently large. Then the argument of w is roughly if 
(actually, ±ip but let us focus on w with positive imaginary part). Then, for R > Rq large 
enough and d < do small enough (though positive) it follows from basic geometry that the 
argument of w{q'^ — 1) can be bounded in {n/2 + b, 37r/2 — b). In order that q^w avoid Ca,,p- 
it suffices that d be less than a constant times \w\~^ and that \w\q^ be small (but still of 
order 1). This implies that R can be chosen roughly as — lug \ w\. 

For the moment we will only make use of the contour Ca,ip- 
Lemma 4.10. Fix k > 1 and consider positive di, . . . , d^- Then, for any ip G (0, 7r/4), 

(-l)^g^ f f T-f za-Zb f{zi)---f{zk) ^ 
^^^ = T^T^k — /■■■/ 11 dzi---dzk, 

with f{z) = exp((g — 1)72;) YliLi 5^ ^'^^ Zj-contours given by Caj,ip where {aj}jLi are 
positive real numbers such that aj < qaj+i for all 1 < j < N — 1 and < niinj{aj}. 
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Figure 4. (a) w-contour Ca,^- The black dots are di, . . . ,a]\f; (b) s-contour 
which depends on the value of w E Ca,tp', (c) The contour is chosen so 
that q^w sits entirely to the left of Ca,ip as is demonstrated here; (d) The image 

of q'^ for s G Vy^. The value of d and R are also chosen so that the argument 
of w{q^ — 1) is contained in (7r/2 + 6, 37r/2 — 6). This amounts to making sure 
that the angle g'^ — 1 makes with the negative real axis is within {—b, b). 



Proof. The formula in (14 .Sp for /i^ involves closed contours for the Zj which can be chosen to 
be concentric circles centered on the real axis and with their left crossing point with the real 
axis a aj (the za contours containing q times the zb contour for all A < B as well as the poles 
at the di's). We now proceed to expand these contours to the infinite contours Caj,:^, one at a 
time. We start with Zi. We may freely deform zi to a pie slice shape with radius r 3> 1 made 
up of the portion of Cai,ip with distance from cti less than r, and then a circular arc through 
the real axis: precisely the contour is {ai + e"^t}o<t<r U {ai + e^"^t}o<t<r U {a + e"^r}_^<o-<<^. 
We would like to replace this by the infinite contour Cai,ip- To justify this observe that 
for z with positive real part, f{z) decays exponentially with respect to the real part of z. 
Since G (0, 7r/4), as r goes to infinity, the contributions of the contour integral away from 
the origin along both the pie slice and the limiting contour become negligible, and hence the 
replacement is justified. This procedure can be repeated for Z2 and so on until the contours are 
as claimed in the lemma. Note that as the original contours were positively oriented circles. 
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and the expanded sections were on the left of the circles, the contours are now oriented so as 
to have decreasing imaginary part. □ 

In order to combine the g-moments fik into a Fredholm determinant it is useful to have all 
of the contours coincide. This, however, involves keeping track of the residues which result 
from such a deformation. First, define the g-Pochhammer symbol and g-factorial as 



n-l 



(a; q)n = Yli'^-q' 



{<i\<i)h 



i=0 



Proposition 4.11. Fix k > 1 and consider a meromorphic function f{z) with N poles 
Oi, . . . , Oat with positive real part. Then setting 

fc(fc-i) 



fJ'k 



we have 
IJ'k 



-k ' 



X 



E 

7 

1 - g)^' 



n 



l<A<B<k 



za - qzB zi--- Zk 



dz\ ■ ■ ■ dzk, 



Ah A: 
A=l'"i2™2 



mi!m2! 



det 



Wiq^* — Wj 



l[f{w,)f{qw,)---f{q'^~'w,)dw^ 

i,j=l j=l 



31 



with the Zj-contours given by Ca^^Lp where {ttj}^^^ are positive real numbers such that 
Oij < qaj+i for all 1 < j < k — 1 and ak < minjjaj}, and where the Wj-contours are all 
the same and given by Ca,ip for any a G (0, minj{aj}). The notation X\- k above means that A 
partitions k (i.e., if X = (Ai, A2, . . .) then k = Y2 ^1)1 ^^'^ ^he notation A = i™-i2™2 . . . means 
that i shows up rrii times in the partition A. 

Proof. Word-for-word repetition of that for Proposition 3.2.1 of |16] . 



□ 

4.3.2. A first Fredholm determinant. Before stating our first Fredholm determinant in this 
derivation, let us recall how such determinants are defined. Fix a Hilbert space L^(X, /i) 
where X is a measure space and is a measure on X. When X = C is a simple smooth 
contour in C, we write L'^{C) where fi is understood to be the path measure along C divided 
by 27ri. When X is the product of a discrete set D and a contour C, fi is understood to be 
the product of the counting measure on D and the path measure along C divided by 27ri. 
Let K be an integral operator acting on /(■) G /i) by Kf{x) = fj^K{x,y)f{y)dfi{y). 

K{x, y) is called the kernel of K. One way of defining the Fredholm determinant of K, for 
trace class operators K, is via the Fredholm series 



det(l + i^)i2(x,^) = 1 + ^- / ■■■ / det [K{xi,xj)]'^.^^Y[dfi{x^). (4.6) 

n=l •'^ ^ i=l 

In fact, if an operator K is such that the above right-hand side series is absolutely convergent, 
then we still write it as det(l + K) as short-hand even if K is not trace class. This is 
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sufficient for our purposes as we will deal directly with the expansion. We will often write 
only det(l + K)l2(^x) for det(l + K)L2^^x,^,)■ 

Our ffist Fredholm determinant formula now follows. 

Proposition 4.12. Fix k > 1 and consider positive di, . . . ,aN ■ Then for any (f G (0, '/r/4) 
and any a G (0, minjjaj}), there exists a positive constant C >1 such that for all \(\ < , 

iCa^l.,^ ) =det(l + K)^.(,^„,,~^^) (4.7) 

ISy ' y/oo / MMt=o(ai,...,ajv;p) 

where the kernel K is given by 

^ C7(^^i)/(gt^i)---/(g"^"V) 

with 

N 



f{w) = exp((g - l)-fw) TT - 



am-w 



Recall that the contour C^^ip is given in Definition 

The left-hand-side of (14. 7p is the g-Laplace transform of the random variable q^^ (with re- 
spect to the Cq exponential) — see Section 3.1.1 of [16]. The g-Binomial theorem is analogous 
to Taylor expansion of the exponential. Because q^'^ < 1 it we can justify interchanging the 
summation (of the series expansion resulting from the g-Binomial theorem) and the expec- 
tation and hence we find that the g-Laplace transform can be written as a generating series 
of g-moments. Using the expression for these moments coming from Proposition 14.111 it is 
an easy formal manipulation of terms in a summation to turn this generating series into the 
desired Fredholm determinant. A little extra arguing shows that these formal manipulations 
are numerical equalities and proves the claimed result. The details are given in Section \T7\\ 

4.3.3. A Fredholm determinant suitable for asymptotics. By using a Mellin-Barnes type in- 
tegral representation and analytic continuation we can reduce our Fredholm determinant 
to that of an operator acting on a single contour. The above developments all lead to the 
following result. 

Theorem 4.13. Fix p a Plancherel (see equation (J^) Macdonald nonnegative specialization 
and positive di, . . . ,dN . Then for all ( ^ C \ IR+ 

^ \ =det(l + K^)^,(,~^^) (4.9) 

where Ca,^ as in Definition \4-8\ with any a G (0, mini{ai}) and any ip G (0, 7r/4). The operator 
is defined in terms of its integral kernel 



(Cg^~;g)c 



k^{wy) = ^ [ T{-s)T{l + s){-Crg^,^iqnds 



where 

exp {-fw{q' - 1)) (q'w/dm; q) 



gw,w'[q j = ; ; I I — jz r , l4.iUj 

q'w-w' [w am]q)oo 
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and the contour Vyj is as in Definition \4.8 



The details of the proof of this result are given in Section 17.21 

4.4. Weak convergence to the Whittaker measure. We are now almost prepared to 
relate the above discussion to the Fredholm determinant in Theorem 14.51 The connection 
relies on the following weak convergence of probability measures result. 

Theorem 4.14 (Theorem 4.1.20 of [IS]). Fix N >1, a drift vector a = (ai, . . . , Otv) G , 
and a time parameter r > 0. For positive di,...,dN and 7 > 0, consider a partition 
A = (Ai, . . . , Aat) distributed according to the q-Whittaker measure (i.e., Macdonald measure 
att = 0) MMj=o(^^i5 • • • ! (^N', p) with Plancherel specialization p determined by 7. For a small 
parameter e > let 

q = e-', 7 = re'^, dj = e-™^ A^ = re'^ - {N + 1 - 2j)e-^ Ine + Tje-\ I < j < N. 

Consider the s-indexed measure induced on {Tj}i<j<N- This measure weakly converges, as 
e 0, to the Whittaker measure ^'M(^ai,...,aN,T) on {Tj}i<j<N- 

4.5. Proof of Theorem 14.51 We may now combine the above results to provide a proof of 
Theorem 14.51 We follow approach taken in proving Theorem 4.1.40 of [16] but due to the 
unboundedness of the contours with which we are dealing, there are a number of extra and 
somewhat involved estimates we must make. These are stated as propositions and proved in 
Sections 17.41 and 17.51 

The proof splits into two parts. Step 1: We prove that the left-hand side of equation fl4.9p 
of Theorem 14. 131 converges to ( e~'"^ ) . This relies on combining Theorem 14.141 

(which provides weak convergence of the g- Whittaker measure to the Whittaker measure) 
with Lemma 18.21 and the fact that the g-Laplace transform converges to the usual Laplace 
transform. Step 2: We prove that the Fredholm determinant expression coming from the 
right-hand side of Theorem 14.131 converges to the Fredholm determinant given in the theorem 
we are presently proving (see 14. 2p . 

In accordance with the scalings of Theorem I4.14[ we scale the parameters of Theorem 14.131 

as 

g = e-^ 7 = re-^ = e-^"^ 1 < A; < 

w = q\ ( = -e^e^''\, Xn = re-^ + (N - l)e-Hne + TNe-\ ^^'^^^ 

Step 1: We assume throughout that m G C with positive real part. Rewrite the left-hand side 
of equation f l4.9p in Theorem 14.131 as 

= (eg(a;g))i 



(f'n^N-n\ I \'"'?'^'^9''/MMt=o(ai,...,ajv;p) 
' y^oo / MMt=o(ai,...,a]v;p) 

where 

= (1 - g)~^Cg^^ = -ue-^^ejil - q) 

and 



{{l-q)x; q)c 
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is a g-exponential. Combine this with the fact that eq{x) — )■ uniformly on Re(x) < 
to show that, considered as a function of T/v, eg{xq) — j- e"""^ uniformly for T/v G M. By 
Theorem I4.14[ the measure on Tat (induced from the g-Whittaker measure on Xn) converges 
weakly in distribution as e — )■ to the marginal of the Whittaker measure WM(aj ... a^^.^) on 
the T/v coordinate. Combining this weak convergence with the uniform convergence of eq{xq) 
and Lemma 15^ gives that 



as g —J- 1, or equivalently e — )■ 0. 

Step 2: Recall the kernel in the right-hand side of equation (14. 9 p in Theorem 14.131 It can be 
rewritten as a Fredholm determinant of a kernel with the variables v and v' (recall w = q") 
as follows: 

dei{l + k(^) = det(l + 

Here the space with respect to which this determinant is defined is that of the contour 
specified in Definition 14. 151 below. The kernel is denoted by to denote the dependence on 
u (through C = —e^e^^ u) and e (through q~^). It is given by 

K{v,v') = ^[ h%s)ds, (4.12) 



where 



h^{s) = r(-s)r(i + s) ( ' -4^6^^^^"'"'^ fr J/''-''^^^ (4.13) 



m=l 

where the new term q'" In q came from the Jacobian of changing wtov and where the g-Gamma 
function is defined as 

^'<^'^(^'^ 

The contour on which this kernel acts is the image of the contour Ca,ip under the map 
X I— luqX. There is a fair amount of freedom in specifying this contour, so we will fix a 
particular such pre-image contour. 

Definition 4.15. Let a = 1 + maxaj then we define the contour C^^ as the image of 
q°' + e^'^'M+ under the map x i— ?■ In^x. This contour is illustrated in Figure\^ 

We will assume that K^^ acts on the contour C^ ^^. Note that as e — t- this contour 
converges locally uniformly to a + e^'^^'^-''M+ as can readily be seen by Taylor expanding the 
map X H- lug X. 

It follows from the above observation that the contour on which the kernel is defined 
converges as e — )■ to the contour a Cq,,^ on which the kernel in Theorem 14.51 is defined. 
Let us likewise demonstrate the pointwise convergence of the integrand in the integral (14.121) 
defining kernel K^^ to that of the kernel Ku- 



FREE ENERGY FLUCTUATIONS FOR DIRECTED POLYMERS 



32 



Consider the behavior of each term as g — t- 1 (or equivalently as e — )■ as g = e ^): 

< 



e""'" (t^-^) (4-14) 



^^^^^ (4.15) 



qs+v — qv ^ _|_ s _ ^. 
^qiy - am) T(v - a, 



Tq{v + s-am) T{s + v-am) 
e^''~' exp i^q^iq' - 1)) ^ e^'"+"'/2_ ^4_X7) 

Combining these pointwise hmits together gives the integrand of the kernel given in f l4.2p . 
However, in order to prove convergence of the determinants, or equivalently the Fredholm 
expansion series, one needs more than just this straightforward pointwise convergence. 

There are four things we must do to complete Step 2 and prove convergence of the de- 
terminants. In proving convergence of Fredholm determinants it is convenient to have the 
contour on which the operator acts be fixed with e. 

In Step 2a we deform Q ,^ to a contour Q^^^^ with a portion Ca^^p^^r (of distance < r to the 
origin) which coincides with the limiting contour Cq^^. 

Then in Step 2b we show that for any fixed 77 > 0, by choosing Eq small enough and tq 
large enough, for all e < Eq and r > tq the determinant restricted to L'^{Ca,ip,<r) differs from 
the entire determinant on L^(C^^^) by less than 77. Thus, at an arbitrarily small cost of 77 
we can restrict to a sufficiently large radius on which the contour is independent of e. 

In Step 2c we show that for any 77 > 0, for e small, the Fredholm determinant of 
restricted to L'^{Ca^^^<:r) differs by at most rj from the Fredholm determinant of restricted 
to the same space. 

Finally, Step 2d shows that for tq large enough, for all r > Tq the Fredholm determinant 
of Ku restricted to L'^{Ca^tp^<,r) differs from the Fredholm determinant of on L'^{Ca^^) by 
at most 77. Summing up the steps, we deform the contour, we cut the contour to be finite, 
we take the e — limit, and then we repair the contour to its final form - all at cost Sf] for 
f] arbitrarily small. 

Step 2a: We must define the contour to which we want to deform ^, and then justify this 
deformation as not changing the value of the Fredholm determinant. 

Definition 4.16. Fix (f G (0,7r/4) and r > 0. For a G M, define the finite contour Ca^^p^^r 
to he {a + te^'"^'-^^^ : < t < r}. The maximal imaginary part along Ca,ip,<r is rsin(y9). 
Define the infinite contour C^ ^^^^ to be the union ofCa,y:>,<r with C^ ^^^^^ and C^^^ _^. Here, the 
contour C^, ,^^^^ is the portion of the contour C^ ,^ which has imaginary part exceeding r sin(<y9) 
in absolute value; and the contour Q ,p is composed of the two horizontal line segments 
which join Ca,ip,<r with C^^^^r- These contours are illustrated in Figure\^ 

Now we justify replacing the contour ,^ by C^^^^^- 

Lemma 4.17. For any r > there exists > such that for all e < Eq, 

det(l + i^:)L2(c.^) = det(l + /C)l2{cj,„J- 
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Figure 5. Left: The infinite contour and the hmiting contour Cq,^. Right: 
The infinite contour C^^^^^ (which we deform from Q ,^). 

Proof. The two contours differ only by a finite length modification. We can continuously 
deform between the two contours. We will employ Lemma 18.31 which says that as long as 
the kernel is analytic in a neighborhood of the contour as we continuously deform then the 
Fredholm determinant remains unchanged throughout the deformation. The only things 
which could threaten the analyticity of the kernel are the poles coming from the left-hand 
side terms of fl4.15p and fl4.16p . On account of the condition satisfied by the contour Vqv 
(see Definition 14. 8p . it follows that these poles are avoided. By choosing e small enough, the 
two contours we are deforming between can be made as close as desired. Taking them close 
enough ensures it is possible then to deform between them while avoiding poles of the kernel 
in V oi v' - hence proving the lemma. □ 

Step 2b: We must now show that we can, with small error, restrict our Fredholm determinant 
to acting on the finite, fixed contour Ca,ip,<r- This requires us choosing r > tq for vq large 
enough, and also choosing e < Eq for Eq small enough. 

Proposition 4.18. Fix ip G (0, 7r/4). For any t] > there exists tq > and Eq > such that 
for all r > ro and E < Eq 



det(l + K)mc^^^ - det(l + irOL2(c.„,<o 



< r]. 



The proof of this proposition is fairly technical and is given in Section 17.51 

Step 2c: Having restricted our attention to the finite, unchanging (with e) contour Ca,Lp,<r we 
may now take the limit of Fredholm determinant on the restricted space as e — )■ 0. 

Proposition 4.19. Fix (f G (0,7r/4). For any t] > and any r > there exists Eq > such 
that for all e < eq 

|det(l + i^:)L2(c„„,<o - det(l + Ku)L^ic^^^,,.)\ < V 
where Ku{v, v') is given by the integral ( li.gj) . 

The proof of this proposition is also fairly technical and is given in Section 17.41 



FREE ENERGY FLUCTUATIONS FOR DIRECTED POLYMERS 



34 



Step 2d: Finally, we show that post-asymptotics we can return to the simple infinite contour 
C 

Proposition 4.20. Fix (f G (0,7r/4). For any rj > there exists tq > such that for all 
r > tq 

|det(l + Ku)l\c^^^,^^) - det(l + Ku)l^(^c^^^)\ < V- 

The proof of this proposition is given in Section 17.31 It is a fair amount more straight 
forward than the previous two proofs and hence is given first. 

Having completed the four substeps we may combine Propositions 14.181 14.191 and 14.201 to 
show that for any rj > 0, there exists > such that for all e < Eq, 



det(l + Kdmc,,,) - det(l + i^«)L^(c.„) 



< 3r/ 



where det(l + K^) is as in the right-hand side of equation fl4.9p in Theorem I4.13[ subject to 
the scalings given in (14. lip . Since t] is arbitrary this shows that 

limdet(l + Kc)lHc,,,) = det(l + Ku)lhc^,,). 

The above result completes the proof of Theorem 14.51 modulo proving Proposi- 
tions I4A81 14191 and lOUl 

5. Details in the proof of Theorem 11.31 

As discussed in Section [21 to finish the proof of Theorem 1 1.3 1 we need to show Theorem 12. II 
For K > 0, Definition 11.21 associates the scaling parameters 6'^, and c'^ which appear in 
the statement of this result. The variable k and 6 are dual in the sense that one could 
instead start with some fixed ^ > and then associated scaling parameters Kg, fg and cg. In 
particular, for 9 = 9'^ one recovers fg = and cg = c*^. In the proof it is more natural to 
parameterize everything by 9 instead of k, so we will do it. 

First we prove the convergence to the GUE Tracy- Widom distribution without boundary 
perturbations, since the proof with boundary perturbations is a small modification of it. 



5.1. Proof of Theorem 12.1( a). We first give explicit expansions for some of the functions 
Definition O Let ^i{z) = -^^ 

Then 

oo 

*(^) = -7e + X1 

n=0 

where 7e is the Euler constant. Hence, 



from Definition 11.21 Let ^{z) = ^ln(r(z)) be the Digamma function and fix ^ G (0, oo 

n + 1 n + z 



Kg 



*'(») = Et^V^. (5-1) 



n=0 ^ 



oo ^ 

fg = 9^'{9)-^{9)=jE + J2[ 



n + 29 



n=0 

oo 



71 -h y n 



1/3 



(5.2) 



c. = ,-v^'W2)^/3=(X:(;^) . (5.3) 
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Under the scahng hmit 

we have to show the following: For as in fll.Sp and a contour := Cq,,^, 

lim det(l + Ku)l2(c^) = det(l - KAi)mr,oo)- 

To show this we start with the kernel (11. Sp . replace r(— s)r(l + s) = — tt/ sin(7rs) and then 
perform the change of variable z = s + w to obtain 

^ ^NG{v)-NG{z)^rN^/3{v-z) 

sin(7r(5 — f)) z — v' 

where 

Giz)=\nTiz)-K^ + rz. 

We will show that the leading contribution to the Fredholm determinant comes for v,v' in a 
A^~^/^-neighborhood of 9. Now let us specify the exact choice for the contour Cy as well as 
the contour along which z is integrated. We choos^ 

C,:={e-\y\ + iy,yeR}. 

Cy is a steep descent path (see the footnote in Section [2?T|) for the function Re(G'(f)). The 
path for z is dependent on v, since it has to pass to the left of, or contain the simple poles 
V + l,v + 2, . . ., see Figure [H] (left). Consider the sequence of points S = {Re(t') + 1, Re(t') + 
2, . . .}. There are three possibilities: 

(1) If the sequence S does not contain points in [9,9 + 3cg^A^~^/^], then let £ G No be 
such that Re{v) + £ e [9-1,9] and we set e = c^^iV"^/^. 

(2) If the sequence S contains a point in [9,9 + 2cg^ N~^/^], then let £ G N such that 
Re{v) +ie[9,9 + 2cg^N-^/^] and set e = Sc^^N^^/^ 

(3) If the sequence S contains a point in {9 + 2cg ^A^-^/^ 9 + 3cg ^iV-^/^], then let £ G N 
such that Re{v) + i e {9 - 1 + 2c^^N-^'^, 9-1 + ?>c^^N~^l^] and set e = ^A^-Vs. 

With this choice, the singularity of the sine along the line ^ + e + iM is not present, since the 
poles are at a distance at least Cq^N~^^^ from it. Then, the path for z is given by 

e 

C, ■.= {9 + e + iy,yeR}u\J B^+k, 

k=l 

and By_^k denotes a small circle (radius smaller than 1/2) around v + k and clockwise oriented. 
If £ = then the small circles are simply not present. The idea behind this choice of the 
path Cz is that the ^-contour consists of a fixed line that is (almost) independent of kernel 

^Theorem 11.171 is stated for ip G (0,7r/4) since one uses the quadratic decay (|7.10p to control the hnear 
term in the bound (|7.13p . For ip = tt/A one gets a Hnear decay instead of (|7.10p whose strength depends 
on the parameter a too, it would not strong enough general a. However, in our case, with a = 6, it still 
works, as can be seen from the bound obtained in Proposition 15.21 The proof could also be adapted to any 
other asymptotic direction < (/? < 7r/4 by simply modifying the path away at a distance greather than 
some (arbitrary but fixed with A^) value Rq (one can not employ any angle if G (0,7r/4) right away from the 
critical point since some steep descent properties are then locally not satisfied). 
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Figure 6. Left: Integration paths (dashed) and Cz (the sohd hne plus 
circles at w + 1, . . . , f + £). The small black dots are poles either of the sine or 
of the gamma function. Right: Integration paths after the change of variables 
Cw (dashed) and Cz (the solid line plus circles at w + 1, . . . ,w + tj, with 
p = p{w) e {1,3}. 



arguments, and an additional number of little circles (i.e., poles) as needed. Moreover, the 
leading contribution of the kernel comes only from the cases where £ = (i.e., situation (1)) 
for which e = Cq^N~^/^. 

Also, we do the change of variable 

{v^ v', z} = {^{w), ^{w'), $(z)} with $(^) ■=e + zCg^N-^/\ 

After this change of variable, det(l + Ku)l'2(c^) = det(l + _R'Ar)L2(c^), the path C„ becomes 
(see Figure |6] (right)) 

C^:={-\y\ + iy,yeR} (5.4) 
and the accordingly rescaled kernel 

= i / 

27ri Jc,:=^-HCi) sm{7r{z -w)cg'^N-^/^) z - w' 

where 

G{w) = ln(r(u7)) + few - Kew^/2. 

In Proposition 15.21 we show that for any w,w' G C^, there exists a constant C G (0, oo) 
such that 

\Kn{w,w')\ < Ce^lM-)l 
uniformly for all large enough. Therefore, 

n 
i=l 
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where the factor n"/^ is Hadamard's bound. From this bound, it follows that the Fredholm 
expansion of the determinant, 

det(l + Kn)l2{c^) = X] ~T / ■ ■ / det{KN{wi, Wj))i<jj<n, 

is absolutely integrable and summable. Thus we can by dominated convergence take the 
N ^ oo limit inside the series, i.e., replace by its pointwise limit, 

lim Kn{w,w') = KAi{w,w') ■= — dz- -, (5.5) 

iv^oo zvri Jg-^i/4oo [z — w)[z — w') 



derived in Proposition 15. ![ i.e., we have shown that 

lim det(l + Kn)l^(c^) = det(l + KAi)mc^)- 

The last part is a standard reformulation, which we report in Lemma 18. 6[ see also [65] . This 
ends the proof of Theorem 12.11 

5.2. Pointwise convergence and bounds. The function G satisfies 

oo oo 

G'{9) = G"{9) = 0, G^^\9) = -2J2^-^^ = -2cl G^^\9) = J2j^,, (5.6) 

therefore G has a double critical point at 9. For the steep descent analysis we need to analyze 
the function g{x, y) = Re{G{x + iy)). It hold^ 

Re(ln r(x + iy)) = E (I - ^ In (^^^^^^) ) " 7ex - ^ ln(x^ + y') 
= 22 ( — TT - 2 ^^"^ ^ ""^^ ^ ^ ln(n)l„>i j - 7eX. 

n=0 V"- + / 

Together with (15. ip and ( 15. 2p we get 

g{x,y) = Re(lnr(x + iy)) + fgx - \^Ke[x? - y^) 



n=0 

It follows that 



^ f{n + 29)x-{x^ -y^)/2 1 ,2 2\ wn. 



(5.7) 



n (r v) ■= ^^(^'^^ = V ( ^ + ~ ^ - ^ + ^ ^ 

ax (^ + ^)' {x + n)^ + y^l ^ 



and 

. ^ dg{x,y) / y y \ .p. 

g2{x,y):= — = >^ 77- — ^-1 — , — T^- — ^ • (5-9) 

•^y ^\[9 + n)^ {x + n)^ + y^J 



n=0 



See for example http : //functions .wolf ram. com/06 .11.19. 0001 .01 
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Proposition 5.1. Uniformly for w , w' in a bounded set ofC^, 

lim Kn{w,w') = / dz- r-. -. (5.10) 

Proof. Consider ty, w' in a bounded set of C^,, i.e., tlie original variables v' of order A^~^/^ 
around the critical point 6. For large enough and w bounded, Re((— w)cg ^A^~^/^) G (0, 1), 
and Cz := $"^((3^) = {1 + iy, ?/ G M}. Using (15. 6p we have the expansion 



NG{^w)) = NG{e) - -w^ + 0{w^N-^/'^) 

3 

-NCmz)) = -NG{e) + \z' - fiez^N^^/^ + ^(z^iV^^/s^ 



(5.11) 



l + 0{{z-wfN-^/'')). (5.12) 



with /ig = G(4)(^)ce V24 > and 

IT ceN^^ 

sm{iT{z — w)c'f)^N~^/'^) z — w 

It is also easy to control the u;'-dependence because |;z — > 1. 

Now we divide the integral over z into (a) | Im(2;)| > 6N^^^ and (b) | Im(2;)| < 6N^^^ for 
some (5 > which can be taken as small as desired (but independent of A^). 

(a) Contribution of the integration over \ Im(2;)| > 5N^^^ . We need to estimate 



^ dz- 



27ri Ji+iy,\y\>5m/--^ sin(7r(2; - w)cq ^A^-V3) z - w' 
From (15. lip . (I5.12p . and the fact that w is in a bounded neighborhood of 0, we have 



(5.13) 



^Mj<0{l) / rf^e^^<'(^(*(°»-^(*(i+^J'»). (5.14) 

Setting e = Cg^N~^^^ and doing the change of variable y = yc^^N'^/^ we obtain 

j'OO 

([5J4D < C(ArV3) / ^y^N{a{e+e,y)-9{em 

Js/ce 

The function g{x, y) := Re(G'(x + iy)) is given in (15. 7p . Finally, in Lemma [5.41 we show that 
the path 6'+e+iM is steep descent for the function —G{z) with derivative of — Re(G(5)) going 
to -00 linearly in Im(5). It then follows that fl5J[5|) is of order N^/^e^s{9fi)-Ng{e+efl)^-c^{S)N 
for some positive constant ci{5) ~ 5^ for small 5. But 

Ng{e, 0) - Ng{e + e, 0) = ^ + 0{N-^'^). 

Thus the contribution of the integration over | Im(2;)| > 5N^/^ is 0{e~^'^'^^'>^) for some positive 
constant C2{5) ~ 5^ for small 5. 

(b ) Contribution of the integration over \ Im(2;) | < 5N^^^ . We need to determine the 
asymptotics of 

_C-liV-l/3 /■ ^g7VG(*(«;))-JVG(<&(2)) ^r{w-z) 

dz — T — — TT^r 7. (5.16) 



27ri yi+ij/,|s/|<5ivi/3 sin(7r(2; - w)cq ^iV-i/3) z-w'^ 
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Using the expansion (15.1 ip and fl5.12p we get 

27ri A_i5jvi/3 (z - w){z - w') e^'Vs-™' 

(5.17) 

Denoting 2; = 1 + iy we have 

Re(zV3) = -3y^ + 1, Re(z^) = / - 6^^ + 1. 

The convergence of the integral is controlled by e"^"^*^ ^^^-3y^ ^ employs the bound 

|e^ - 1| < Ixlel^^l with x = Oiw^N-^/'"^; z^N"'^/'^) to control the error terms. Altogether they 
are only of order 0{N-^/^), i.e., we have obtained 

flETTD = OiN~^'^) + — / dz-- — --^7^ 

Finally, we deform the integration contour to the following one: from 5N^'^{1 - i) to 
SN^/^{1 The error term is again of order e"'^^^^)^. However, with the new contour, 

using again |e^ — 1| < |x|el^l but with x = —figz^N~^/^ one sees that the eliminating the 
quartic power in z amounts in an error of order 0{N^^^^). The last step is to replace S by 
00 in the integration boundaries. This leads to an extra error O^e"^^^^"^^) with some positive 
constant 03(6) ~ 6^ for small S. 

To summarize, we first choose S small enough so that all the Ck{S) > 0. Then for all 
large enough we have shown that the contribution of the integration over | Im(z)| > SN^^^ is 
of order 0{e~^^^^^^) and the integration over | Im(2;)| < 5N^^^ is given by 

_i pe^'/'^co 1 ^z^/3-rz 

0{N-^/^) + — / dz- 



2ni Je-^^/^oo {z -w){z- w') e"''/3-ru, • 
Taking the N ^ 00 limit we obtain the result. □ 
Proposition 5.2. For any w,w' G C^, there exists a constant C G (0, 00) such that 

\Kn{w,w')\ < Ce^lM-)l 

uniformly for all N large enough. 

Proof. Since the z-contour can be chosen such that \z — w'\ > 1/2, we can estimate the 
absolute value of the factor [z — w')~^ by 2 and discard it from further considerations. For w 
in a bounded set of Cw, the statement is a consequence of the computations in the proof of 
Proposition 15.11 Thus, it is enough to consider w = —\y \ + iy for y > L, for L which will be 
chosen large enough (but independent of A^). In the original variables v,v', this means that 
we need to consider v = 9 — \y \ + iy for y > Lcg^N~^^^ . Let v = $(w), v' = $(w'), then the 
kernel Kn is given by 

^N{G{v)-G{e))+r{v-e)ceN^^'^ r ^^^NG(e)-NG(z) ^r{e-z)cgN^/^ 

Kn(w,w') = — 7^ ; / dz 



CqN^/^ 2rc\ Jq. sin(7r(2; — v)){z — v') 

We divide the bound dividing in two contributions: (a) integration over 6 -\- e + iM, with 
e = pCg^N~^^'^ (with p G {1, 3} depending on the value of v, see the proof of Theorem 12.1( a) 
above), and (b) integration over the circles B{v + k), k = 1, . . . , i{v). 
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(a) Integration over 6 + e + iM. The relevant dependence on v is in the prefactor 
^N{G{v)-G{e))+r{v-e)c0 sine. The dependence of e on v is marginal, as the needed 

bounds can be made for any e small enough. The estimates as in the proof of Proposition 15. II 
imply that this contribution is bounded by 

(j^N{Re{G{v))-G{e))+r{Rc{vye)ceN^^'' ^ ^N{g{e-y,yyg{e,0))-rycgN^/^ ^ 

where we used the parametrization v = 9 — \y \ + iy and, by symmetry, considered only y > 0. 
In Lemma [5. 3 1 we show that g{6 — y, y) is strictly decreasing as y increases and for large y the 
derivative goes to — oo (logarithmically). Thus, for any fixed 6 > 0, there exists a constant 
Ci > such that for all y > S, dyg{6 — y,y) < —Ci. In Lemma [5. 3 1 we also show that for small 
y, g{6 — y,y) — g{9, 0) = —^c^y^ + 0{y'^). Therefore, we can choose 6 > small enough such 
that: 



3 



(1) for Lcg <y<S,g{9- y, y) < g{9, 0) - ^c^y 

(2) for y > 6/2, dyg{e -y,y)< -2ci. It follows g{e -y,y)< g{e, 0) - c^y for all y>5. 
Replacing y = lm{v) = lm{w) / (cgN^^^) we get the bounds: 

(1) for L < lm{w) < 5ceN^I^, 

Im(ui)'^/3— r Im(ui) ^ Im(ui)'*/6 ^ g^g— Im(ui) 

for L large enough (depending on r only). 

(2) for Im(w;) > 5ceN^/'\ 

(j^-lra{w){N'^/'''ci/cg+r) ^ ^g- Im(iy) 

for large enough. 

(h) Integration over the circles B{v + k), k = l,...,i{v). This happens only if 
y + Scq'^N^^^^ > 1, where y = Im(f) = lm(w)/ (ceN^^^). The contribution of the integra- 
tion over -B^+fc (up to a ± sign depending on k) is 

^NG{v)-NG(v+k) ^-rkcgN'^/^ 

{v + k — v') 

We have \v + k — v'\ > l/\/2, thus the contribution from the pole at f + is bounded by 

2^N{g{e-v,v)-g{e-v+k,v))^\r\kceN^/^ _ 

Define the function h{v, k) := g{d — v,v) — g{6 — v + k,v). In Lemma 1531 we show that h{v, k) 
is strictly decreasing as a function of k, for G [0,i/ + e] (we have a positive 6 instead of e, 
but for large enough e < 6). Also, k < i{v) = [y + e\, so that the contribution of the 
poles at V + 1, . . . ,v + £{v) is bounded by 

2£(t;)e^'^("'^)+l"l^('')"«^'^'. (5.18) 

We consider separately the cases {!) y < (i.e., Re{v) > 0) and (2) y > (i.e., Re{v) < 0). 

(1) For y < 6, from the bound on dkh{v, k), see Lemma [5.5[ we get h{v, 1) < —ycl/4. 

(2) For y>6,we know that h{v, 1) < for all y and when ?/ — )■ oo, dkh{v, k)\k=o — —yi^e- 
Since the function h{y, 1) is continuous in y, there exists a positive constant c > 
such that h{y, 1) < —cy for all y > 9. 
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Thus, with c' = min{c, Cg/4} we get 

for large enough. This ends the proof of the Proposition. □ 

Finally let us collect the lemmas on the steep descent properties used in the propositions 
above. 

Lemma 5.3. The function g{6 — y, y) is strictly decreasing for y > 0. 

For y ^ oo it holds dyg{6 — y,y) — hi{y). 

Fory\0 we have g{9 -y,y) = g{9, 0) - fc^y^ + 0{y^) . 

Proof. Using (15. 8p and (15.91) we have 

dg{0-y,y) \ fa \ 

- = 92{0-y,y) -gi{e-y,y) 



dy 

(X 

El _L u-r u - ^y \ _ _ \ - 



1 n + e-2y \ ^ 2y^ ^^'^^^ 



^e + n {n + e- yf + y^ J {0 + n){{e + n - yf + y^) 

which is for y = and strictly negative for y > 0. The asymptotics for large y are obtained 
by writing (I5.19P as 

U>{-y + e + iy)- U>{-y + e-iy)- ]^^{-y + e-iy)- ]^^{-y + e + iy) + ^>{e) 
and using the large-z expansion 

^{z)=\Yi{z)- — + 0{z-'^). (5.20) 
Taylor expansion gives the small y estimate. □ 

Lemma 5.4. For any x > 6, the function g{x, y) is strictly increasing for y > 0. 
For y ^ oo it holds dyg{x, y) ~ Kgy. 

Proof. From (15. 9p we have 



9y \(^ + '^)^ {x + ny + y'^ 

which is for y = and for y > is strictly positive. For large y, the second term goes to 
zero, leading to the estimate. □ 

Lemma 5.5. Let y > be fixed. The function 

Ky^ k) ■= 9{0 - y,y) - g{9 - y + k, y) 

satisfies h{y,0) = 0, h{y,k) is strictly decreasing for k G [0,y]. 

For any S G (0, 9), y > S, h{y, k) is strictly decreasing in k E [0, ?/ + 5/2]. 

Fory -> oo, dkh{y,k)\k=o ~ -yf^e- 

Fory< 9, dkh{y,k) < - '^ for k E [0,y]. 
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Proof. From flS.Sp we have 

dh{y, k) f9 + n + y — k 6 + n — y + k 



-gi{e-y + k,y) = " ( 

n=0 ^ 



_ ^ (^ + n)(y^ - (fc - + {y- kf + 1/2(1/ - k) 



re=0 



which strictly negative for k G [0,?/]. 
The second statement follows from 

{e + n){y^ -{k- yf) + {y - kf + y\y - k) 

>e{y^-{k-yf) + {y-kf + y\y-k) 

> 6(6^2 - 6^4) + V/2 - 5^/8 - 6yy2 > 6^^ 

To get the asymptotics of the derivative for large y, we can rewrite 



dh{y, k) 



dk 



= ^{6) - ^\e)y - \^{-y + e-iy)- ]-^{-y + e + iy) 



and use fl5:20D and *'(^) = kq. 

Moreover, for k G [0, y] and y < 6 we have the bound 

dh{y,k) {y^-{k-yf) 



9k ~ ^^{e + n){{e + n-y + ky + y^) 



- ^^^^ie + n)ii9 + n-y + kr + y^)- ^^Sj2(^ + n) 



ky4 



(5.21) 



□ 



5.3. Proof of Theorem 12.1( b). Now we turn to the proof of the Theorem with boundary 
perturbations. 

Note that due to the ordering of the a^'s, 6i > 62 > ■ ■ ■ • Call b = hi. The scaling of the Oj's 
implies that the contours and can be chosen as before except for a modification in a 
A^~^/^-neighborhood of the critical point, since they have to pass on the right of 6 + hc^^N~^/^ 
(see Figure [7]) . Let us denote 

Pin^ z a) - r(^H-^) r($(.)) 

P[w,z,a).- ^^^^^^^ r($(^)-a)- ^^-^^^ 

Then, the only difference with respect to the kernel (]5.10p is that in the -t- 00 limit there 
might remains a factor coming from YYk=i ^k)- 

Using ~ z'^-^l + 0{l/z)) (see (6.1.47) of [1]), for any w,z on C^,Cz we have the 

bound \P{w,z,ak)\ < Ce'=l'^fel(|i"^(«')l+|ii°WI)^"'/' for some constants C, c. The local modifica- 
tion of the paths has no influence on any of the bounds for large w and z, so that the proof 
of pointwise convergence and of the bounds are minor modifications of Proposition 15.11 and 
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^ C 



• • • • ^« (?) (5 




0{1) 

Figure 7. Perturbation of the integration paths, compare with Figure [6] 
(right). The white dots on the right are the values of bi, ... ,bm- 



Proposition I5.2[ It remains to determine the pointwise hmits of P{w, z, at) as N oo. 
Case 1: If limsup^_^oo('^fc(^) ~ 0)N^/^ = — oo, then 

lim P{w, z, ttk) = 1. 

Case 2: If limsup^^^(afc(iV) - e)N^/^ = b^, then 

r(<l>H - $(6,)) r($(z)) 



lim 



N^oo r($(w)) r($(z) - <i)(6fc)) 

T{{w - bk)c^^N-'/^) T{e + zc^^N-'/^) 



lim 



z -bi 



N- 



oo T{e + wCs^N'y^) T{{z-bk)c-^N-y^) w-bk 



because T{z) = z ^ — ■Je + ^{z) as 2 — )■ 0. 
Therefore, one obtains 



where 



lim det(l + Kn)l-^(c^) = det(l + i^BBP,b)L2(c„ 



n 



i^BBP,fe(s, s') = — / dz 

27riJe-W4oo [z - w){z - w') ^J-^w - bk 

The reformulation of Lemma 18.71 complete the proof. 



(5.23) 



6. Details in the proof of Theorem 11.101 

As discussed at the beginning of Section [5l in the proof we parameterize using the position 
of the critical point 6 instead of n. Let us set T G (0, oo) and consider the scaling limit 

e := y/N/T, T = KgN, u = Se-^fo. 
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One has the following large-^ expansion of (15. ip and ( 15. 2 p gives 

A = 1 - !n(») + l + ^ + 0{e-% 

Thus, 

T = KgN = y/TN +\T+0{N-^/^), 

Equivalently, we can set r = VTN, then 6 = ^/N/T + | - ^^/T/N + 0{N~^/^), so that 

T = T- T^/^N^/^ + i|TViV + C(iV"3/2), 

^ ^ g^~N^^NlniT/N)~lVTN+T/i\+OiN-^) 

As shown in Section |3l what it remains is to prove Theorem 13.31 We first prove the 
statement for the unperturbed case, and then we will show how the generalization is obtained. 

6.1. Proof of Theorem 13.3( a) . We have to determine is limAr^oo det(l + -f^u)L2(c«)- Con- 
sider the case of a drift vector 6 = 0. The path C^, is chosen as 

C^ = {e- 1/4 + ir, |r| < r*}U{9e'\t* < \t\ < n/2}U{9- \y\ +iy,\y\ > 9}, 

where r* = y/9/2 - 1/16, t* = arcsin( ^^1/2^ - 1/16^2). The path is set as 

e 

C, = {9 + p/A + iy,y eR}u[jB,+k, (6.1) 

k=l 

where Bz is a small circle around z clockwise oriented and p G {1,2} depending on the 
value of V, see Figure El More precisely, for given v, we consider the sequence of points 
S = {Re{v) + 1, Re(t') + 2, . . .} and we choose p = p{v) and £ = l{v) as follows: 

(1) If the sequence S does not contain points in [6', 6* + 1/2], then let ^ G No be such that 
Re(u) + £ G [9-1,9] and we set p=l. 

(2) If the sequence S contains a point in [6*, 6' + 3/8], then let £ G N such that Re(f ) + £ G 
[^,^ + 3/8] and set p = 2. 

(3) If the sequence S contains a point in [9 + 3/8,^ + 1/2], then let £ G N such that 
Re(w) +Ee[9- 5/8, 9 - 1/2] and set p = 1. 

With this choice, the singularity of the sine along the line ^ +p/4 + iM is not present, since 
the poles are at a distance at least 1/8 from it. Also, the leading contribution of the kernel 
will come from situation (1) with £ = and p = 1. We denote 

a := (2/r)^/=^ 

and we do the change of variable 

{v,v\z} = {^{w),^{w'),^{z)} with ^{z):=9 + za 
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Figure 8. Left: Integration paths C„ (dashed) and (the sohd hne plus 
circles at w + 1, . . . , w + where ^+ = ^+p/4 and 6'„ = 6' — 1/4, the small black 
dots are poles either of the sine or of the gamma function. Right: Integration 
paths after the change of variables Cw (dashed) and Cz (the solid line plus circles 
at w + 1, . . . , w + £), with -p = p{w) G {1,2} 



and 

Ke{w, w') := aK^^w), ^w')) = — / dz-—. . (6.2) 

27ri Jq_ sm(7r(2; — w)a) z — w' 

After this change of variable, the paths = $~^(C,t,) and Cz = $^^(Ci) are given by 

C ={-1/(4^7) + ir/a, \r\<r*}U {(e" - l)9/a, t* < \t\ < n/2} U {-\y\ + iy, \y\ > O/a}, 

e 

Cz ={p/(4a) + iy, 2/ G M} U B^u^k/a, 

k=l 

(6.3) 

where r* = ^9/2 - 1/16, t* = arcsm 

(^1/2^-1/16^2)^ and 5, is a small circle around z 
clockwise oriented. After this change of variable, we have 

det(l + K^)lHc.) = det(l + Kg)L2^c^y (6.4) 

Thus, we need to prove that 

lim det(l + Kg)L2{Cu,) = det(l - Acdrp)l2(r+) 

with A'cDRP given in Definition II. Ill The proof is very similar to the second part of the proof 
of Theorem 12.1( a). where this time the convergence of the kernel is in Proposition 16.11 and 
the exponential bound in Proposition 16.21 We then obtain 

lim det(l + Ke)L2(c^) = det(l + A^cdrp )l2(c») 

with A'cDRP given in (16. 6p . Lemma 18.81 shows that the limiting Fredholm determinant is 
equivalent to det(l — A'cdrp)l2(r+) and thus completes the proof of Theorem 13.3( a). 
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6.2. Pointwise convergence and bounds. The leading contribution of the Fredholm de- 
terminant and in the kernel comes from w, w', z of order 1 away from Q ~ (9(-\/iV). The scale 
for steep descent analysis is A^^ instead of as in the case of the convergence to the GUE 
Tracy- Widom distribution function. So, the function whose real part has to be controlled is 
this time 

~ _ G{e + ez) 
C^[Z) .- — - — , 

that satisfies ^ 

G(3)(o) = -1 + C)(r^), 

G(^)(0) = 2 + O(r^), 

~< , (6.5) 
G(")(0) = C(l), n > 3, 

G'("^(^) = r"+^G(")(0). 

For asymptotic analysis we need to control the real part of G, which we denote 

g{e + ex,eY) 



g{XX) :=Re(G(X + iF)) 



e 



In Lemmas 16.51 16.41 and 16.61 we will analyze the steep descent properties for G (those are 
analogs of Lemmas 15 . 3tl5 . 5 p . that we use to prove Proposition 16. II and Proposition 16.21 below. 

Proposition 6.1. Uniformly for w,w' in a hounded set of 0^,, 

lim Kg{w,w') = KcDRp{w,w') 

where 

~ —If (TttS'^^"'")'^ ^z'-^/a-w^/a 

KcDRpiw,w') = — dz —. (6.6) 

27ri J j_ sm(7r(z — w)a) z — w 

Proof. First remark that the only dependence on in the kernel (16. 2p is in the factor 

exp [N{G{<l>{w)) - Gmz)))] = exp InO (G{wa/e) - G{za/e) 



Let w, w' be in a bounded set of around the origin. For N large enough and w bounded 
in C^, Re{wa + 1) > 1/2 and Re{{z — w)a) G (0, 1) so that we have £ = and p = I, i.e., in 
this case Cz = {j^ + iy , y E M.} . We have 

7VG'($(w)) = NG{9 + wa) = NG{9) + —G^^\9)a^w'^ + 0{Nw^/d^) 

6 

= NGi9) - + 0{Nwye\ Nwye^) (6.7) 

= NG{e) + o{wye) 

3 

where the ^-dependence in the error term follows from G^^\6) = 0{9^^) and then we used 
the expansion (16. 5 p for G^^\9). 

We divide the integral over z into two parts: (a) | Im(z)| > 9^^^ and (b) | Im(z)| < 9^^^. 
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(a) Contribution of the integration over I Im(^)| > ^^/^ 

w,w' on of order 1 and 

z e C„ \z - w'\ > C(l), I sin(7r(^ - w)a)-^\ = So, 

dYexp [Ne {g{0,0)-g{{Aaer\Y))] . 

-2/3 

From Lemma [6.41 we have that —g{{4a6)~^, Y) is strictly decreasing with derivative going to 
— oo as y goes to infinity. Then, the integral over Y is bounded and of leading order 

exp [Ne {giO, 0) - m^0)~\ r^/^))] . (6.8) 
The estimates for small Y of Lemma [6.41 with X = (4o"6')~^ and Y = 9"'^^^ lead then to 



dnH) < exp 



N9 (giO, 0) - gii^a9)-\ 0) - ^ + 0(r ^^3; 



C(l)exp 



12 



where we also used the fact that g{{4:a9) ^,0) = g{0,0) + 0{9 Thus, the contribution 
in the kernel from the integration over | Im(2;)| > 9^^^ is of order e'^^^^'^ for some positive 
constant c > 0. 

(b) Contribution of the integration over \ Im(z)| < 9^^^. We need to estimate 



27ri 



1 f (T7rS'(^~"')'^ ^NG{<l>{w))-NG{^{z)) 

^+iy,\y\<ey3 sm(7r(2; - w)a) z - w' 



Unlike the scaling where we have proven the convergence to the GUE Tracy- Widom distri- 
bution, in this case the sine function survives in the limiting expression and we do not have 
to employ the quartic term in the estimates (since it was used only to control the error term 
of the sine). 

First we verify that the convergence is controlled by the third order term. For this purpose, 
we set z = iy + l/(4cr). Then, using ( 16. 5 p we obtain (as in ( 16. 7p ) 

-NG{9 + az) = -NG{9) + — + 0{z^l9). 

3 

The real part of the cubic term is given by 

In our situation we have \y\ < 9^^^, therefore — ^ dominates 0{z^/9) for large 9 (since 
We have 

(EH) = — dz- 



2n\ 



We divide the integration in (b.l) 9^/^ < \y\ < 9^^^ and (b.2) \y\ < 9^1^. Since the qua- 
dratic term in y from (16.101) dominates the others, the contribution of (b.l) is only of order 
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0{e '=1^^'''') = 0{e '^2^^'"') for some constants Ci,C2 > 0. The contribution from (b.2) is given 
by 

1 r (T7rS'(^~"')'^ ^z'^/s-w-^/s+oiw^/e^zye) 

dz^— . (6.11) 



27ri 7j_+ij^,|y|<ei/6 sm(7T(z-w)a) z - w' 



For \y\ < e^/^, 0{zye) = Using |e^' - 1| < |a;|el^'l for x = 0{z^/e) and then for 



X 



Thus, 



0{w'^/6) we can delete the error term by making an error of order _ Q^p^^ 



1/6 ^ 



J ^+iy,\y\<e^/'' 



Finally, extending the last integral to ^ + iM we make an error of order 0{e~^^^ ) for some 
constant C3 > 0. 

Putting all the above estimates together we obtain that, for w, w' G Ci„ in a bounded set 
around 0, 



Ke{w,w') = 0{N"^/^) + -^ [ dz 



_1 f anS^^^^''"' e^V3--i«'^/3 

sin(7r(2; — w)a) z — w' 



□ 



Proposition 6.2. For any w, w' in 0^, uniformly for all N large enough, 

\Ke{w,w')\ < Ce"l^'"("')l 

for some constant C . 

Proof. First recall the expression of the kernel, 
Kf)(w,w') = — — / dz- 



2-Ki Jq_ sin(7r(2; — w)a) z — w' 

(^-wa^NGmw))~NG{e) ^ I ^T^^ e 



2ni In sin(7r(2; — 1x1)0] 



As in the proof of Proposition 15.21 the dependence on w' is marginal because (a) we can 
choose the integration variable z such that I2; — > 1/(4(t) and (b) we will get the bound 
through evaluating the absolute value of the integrand of (16. 2p . 

Case 1: w G {-1/(4(t) + iy,\y\ < r*/cr} with r* = ^/9/2 - 1/16. In this case, the 
integration path for 2; is 1/ {4a) + iM and no extra contributions from poles of the sine are 
present. The factor 1/ sin(7r(2; — w)a) is uniformly bounded from above. Doing the change 
of variable 2; = + i— we get 



\Ke{w,w')\ < 0(i)e^i^'=(f?(^+--))-^f?W / rfye^^(^(°'°)-^(^"'^))^ (6.12) 

Jr 

with e = 1/ (4:6). The estimates as in the proof of Proposition 16. II on the integral over Y yield 
( 1632|) < C(l) X e^^<c;is+^'^))-NG{e) ^ ^ ^NeRc{G{wa/e))-NeGio) _ (^q -^^-^ 
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Since \wa/6\ < 0{9 ^/^) is small, we can use Taylor expansion and with (16 .Sp we obtain 

NeG{wa/e)) - Nedio) = -Iw^i + o(ri)) + ^w^'ii + o{e-')), 

substituting w = — l/(4cr) + iy and taking the real part we get 

Ne Re{G{wa/e)) - NeG{0) = -^y' + + 0{l) + 0{y'/e, y'/d'). 

Now, for \y\ < -^'Oj^jo, ^y^ < j^y'^ and the quadratic term dominates 0{y^ / 9,y^ /O"^) for 
large 9. Therefore, for all 9 large enough, we have 

N9 Re{G{wa/9)) - N9G{0) < -^y^ + C(l). 

OO" 

Consequently, 

\Ke{w,w')\ < 0(l)e-8^l^"^("')l' < Ce-I^'"^'")! 
for some finite constant G. 

Case 2: w G {(e'* - l)9/a,t* < \t\ < 7r/2} U {-\y\ + iy,\y\ > 9/a}. We divide the 
estimation of the bound by dividing into the contributions from (a) integration over ^ + iM 
with p G {1,2} depending on w (see the definitions after (16. ip ) and (b) integration over the 
circles Byj+k/a, k = 1,. . 

Case 2(a). First notice that the estimate (I6.12p of Case 1 still holds with the minor 
difference that e = p/{4:9) where p G {1,2} depending on the value of w. Then, also (I6.13P 
still holds, so that we need only to estimate N9'Re{G{wa/9)) — N9G{0). 

For w G {(e'* — l)9/a,t* < \t\ < 7r/2}, in Lemma [6.31 we show that 5f(cos(t) — l,sin(t)) — 
5(0,0) < — sin(t)''/16. Replacing lm{w) = sm{t)9/a and using |Im(u;)| > ^^6*72 — 1/16 we 
obtain 

N9Re{G{wa/9)) - N9G{0) < -Ci\lm{w)\y9 < -C2I Im(w;)|v^ < -|Im(w)| 

for all 9 large enough, where Ci,C2 are some (explicit) constants. This is the desired bound. 

For w G { — |y|+iy, \y\ > 9/a}, from Lemma 1631 it follows that there exists a constant C3 > 
such that (9y5'(— y, Y) < — C3 and from Lemma [6.31 we know that g{—l, 1) —g{0, 0) < —1/16. 
Thus, for C4 = min{cr/16, C3} it holds 5f(— 1, 1) —^(0,0) < —C4Y for all \Y\ > 1/a. This 
means that 

N9Re{G{wa/9)) - N9G{0) < -c^N9\ lm{w)\/9 < -\ lm{w)\ 

for large enough, giving us the needed bound. 

Case 2(b). It remains to check that the extra contributions of the poles of the sine also 
tend to zero exponentially in | Im(w)|. The contribution of the integration over Byj+k/a is (up 
to a ± sign depending on k) given by 

igk^NG{i5>{w))-NG{^iw+k/a)) 

w + k/a — w' 

Let us set h{Y, k) := g{—Y, Y) — 'g{—Y + k, Y). From Lemma [6.61 it follows that the largest 
contribution comes from the integration over By^jj^ij^j. We have at most 0(|lm(u;)|) poles 
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and also \w + k/a — w'\ > 0{l/6) (the worst case is at the junction between the arc of circle 
and the straight lines). Thus, the contribution of all the poles is bounded by 

0{e\ Im(w;)|)5li'"('")le^^(*("))-^^(*("'+i/'^)) = 0{9\ Im(u;)|)5l^'"('")le^^'^(^'i), (6.14) 
where Y = \ lm(w)\a/6. We consider separately the cases l/\/26 <Y<1 and Y > 1: 

(1) For 1/V29 < F < 1, the bound on dkh{Y, k) leads to h{Y, 1) < -F/8. 

(2) For r > 1 we^know that h{Y,l) < and dkh{Y,k)\k=o ~ -Y as Y^-^ oo. By the 
continuity of h{Y, 1) in F, there exists a constant C5 > such that h{Y, 1) < —c^Y 
for all Y>1. 

Therefore, with cg = minjcs, 1/8} and inserting Y = \ lm{w)\a/6 we have 

(Km> < C(e|Im(^/;)|)^l^"^('")le"^"^"l^'"("')l < 0(l)e"l^"^('")l 

for large enough. We have shown that also the contributions of the poles have the desired 
bound. □ 

Lemma 6.3. The function 'g{cos{t) — l,sin(t)) is zero at t = and strictly decreasing for 
t e (0,7r/2]. 

Fort G [0,7r/2] and 9 large enough, dtg{cos{t) — l,sin(t)) < — sin(t)(l — cos(t))/2 so that 

g{cos{t) - 1, sin(t)) - 5^(0, 0) < - sin(t)Vl6. 
Proof. We have g{cos{t) — l,sin(t)) = 6~^g{6 cos{t),6 sm(t)), thus 

gfi^(cos(t) - l,sin(t)) ^ cos(t)(72(gcos(t),gsin(t)) -sin(t)(7i(gcos(t),^sin(t)) 
ot 

^ 2g2sin(t)(l -cos(t))(2r2cos(t) + g) ^^"^^^ 

~ {2ne cos{t) + + e^){n + ey 

is strictly negative for all t G (0, 7r/2], which shows the first result. We can further bound 

26^ 



( Km < - sin(t)(l - cos(t)) 



n=0 



{2ne cos{t) + n"^ + e^){n + e) 



2 



20'^ 

< -sin(t)(l - cos(t)) ^ = - sin(t)(l - cos(t))(| + 0(6-')). 

Thus, for large enough Q, the derivative is bounded by — sin(t)(l — cos(t))/2, t G [0,7r/2]. 
Integrating over [0, t\ gives 

5^(cos(t) - 1, sin(t)) - 5^(0, 0) < -(1 - cos(t))V4 < - sin(t)Vl6 
for t G [0,7r/2]. □ 

Lemma 6.4. For any X > 0, the function g{X, Y) is strictly increasing for Y > 0, with 

dYg{x,Y)>dYmy)- 

For Y \ 0, dygiO, y) = + 0{Y^/e- Y^), so that 

giX, Y) > giX, 0) + Y^U + OiY^e; Y'). 
For Y -^00 It holds dygiX, F) ~ F . 
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Proof. We have 



OO / 

g^{e^ex,eY) = Y,[-^ 



dY ' ' ^{(e + nV (6 + ex + n)^ + e'^Y'^ 

which is for F = and for F > is strictly positive. The inequahty 



dyg{X,Y)>dyg{^,Y) = Y, 



whose expansion for small Y and large Q is given by F'^/S + 0{Y'^ jQ] Y'^\ 

For large F, the second term becomes irrelevant with respect to the first, so that 
dYg{X,Y)^eKQY = Y{\\0{B-^y). □ 

Lemma 6.5. The junction g{—Y, Y) is strictly decreasing for Y > 0. 

For F ^ oo It holds dygi-Y, F) ~ - ln(F). 

ForY\0 we have g{-Y, Y) = g{0, 0) - ^Y^ + 0{Yye; Y^) . 

Proof. As in the proof of Lemma 15.3^ we use flS.Sp and fl5.9l) to obtain 



dY ^.v^ (^ + ^)((^ + ^_^y)2 + ^2y2)• 

which is for Y = and strictly negative for F > 0. 

We can rewrite the sum with the variable rj = n/6. Then, for large (but still fixed) 6 
the sum over r] G {0, 1,2, . . .}/6 is very close to the integral over t] G [0, oo). From this one 
deduces that for large F, dygi—Y, F) ~ — ln(F). The asymptotics for F \ is obtained by 
writing the Taylor series of G{Z) around Z = 1 and taking the real part of it. □ 

Lemma 6.6. Let Y > be fixed. The function 

h{Y, k) := g{-Y, F) - g{~Y + k, F) 

satisfies h(Y, 0) = 0, h{Y, k) is strictly decreasing for k G [0, F] . 

For any 5 G (0, 1), Y > 5, h(Y, k) is strictly decreasing in k & [0, F + 5/2]. 

ForY^oo,dJi{Y,k)\k=o---Y. 

For Y < 1, dkh{Y, k) < -kY/A for k G [0, F]. 

Proof. The first statements follows directly from Lemma 15.51 The asymptotics for large F 
can be obtained by approximating the sums in dkh by integrals. The bound for F < 1 follows 
from fICT]) with k ^ 9k, y ^ 9Y, and J2n=oi^ + ^)"^ > l/{29^). □ 

6.3. Proof of Theorem 13.31 (b). Now we consider the perturbed case, where 

ak:=9 + bk, k = l,...,m. 
Then, the change of variable as in (16. 2p leads to the kernel 

Ke{w,w') ■.= aK^{^{w),^{w')) = — / dz-— — l[P{w,z,b,) 

/vri Jq^ sm(7r(z — w)a) z — w 



9Y, 9Y) -gi{9- 9Y, ^F) = - ^ 
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0{1) 

Figure 9. Perturbation of the integration paths, compare with Figure M 
(right). The white dots on the right are the values of bi/a, . . . , bm/o-. 



where the perturbation term is 



P{w,z,k 



•k 



The difference from Theorem 13.31 (a) is that now (as it was the case for Theorem 12.1( b)). the 
paths Cz and C^, have to be locally modified around the critical point, 6, so that they remains 
on the right of all the bi/a, . . . , b^/a, see Figure [9] for an illustration. 
We just have to show that 

lim det(l + Ke)L^{Cn,) = det(l - i^cDRP,b)L2(R+) 



with KcDRP,b as in (11. 6p . 

The proof is a minor modification of the one of Theorem 12.1( b). The local modification 
of the paths have no influence on the bounds for large z and/or for large w. This is because 
NG{9 + bk) — NG{6) = 0{1) and the path for z is the same away from a distance 0{1) from 
the origin. What remains to be clarified is the limit kernel. We can choose the path to 
be as before with a small perturbation (e.g. a circle) around so that it passes on the right 
of the all the b^s. Then, we modify the path in the same way, i.e., by doing the same 
small perturbation but shifted to the right by l/2cr (see Figure IH] too) . This ensures that we 
do not get extra poles from the sine. Finally, for the pointwise convergence of the kernel the 
new term remaining is 

lim P{w,z,bk) = ry- 

A^^oo 1 [az — bk) 

Indeed, since $(z) = 6 + 0{1) and ^{w) = 6 + C(l), then lim^^^o r{|{g)gTI = 1- 

Finally one reformulate the Fredholm determinant into one on L^(M_|_) in the same way 
as the unperturbed case of Lemma 18.81 The only small difference is that we the first step 
requires Re(z — w') > 0, which holds only for bk < j, for all k. Under this condition the 
rewriting holds. By looking at the final expressions one verifies that both sides are analytic 
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in the parameters 61, ... , km- Thus we have equahty by analytic continuation. This ends the 
proof of Theorem 13.3( b). 

7. Details in the proof of Theorem 14.51 

7.1. Proof of Proposition 14.121 This closely follows the proof of [16] Proposition 3.2.8 
and Corollary 3.2.10. However, in that case the contour playing the role of Cq, ,^ is bounded 
whereas it is unbounded presently. As such, some additional estimates must be made, so we 
include the entire proof here. 

First observe that we may combine the g-moments /i^ = {q^^'^)MMt^oia.i,-.-,aN;p) (see defini- 
tions in Section 14. 3p into a generating function 

^qK",) \y /MMt^o(ai,...,ajv;p) 



A:>0 



where kg\ = {q; q)n/{^ — qY ^'^d (a; q)k = (1 — a) ■ ■ ■ (1 — aq^^^) (when k = 00 the product 
is infinite, though convergent since |g| < 1). The convergence of the series defining Gg{() 
follows from the fact that g^'^^ < 1 and 

(C/(i-g))' 



which shows geometric decay for large enough k. This justifies writing 

(C/(i-g))\..A.\ _/ 1 



''=>0 / MM,=o(ai,...,a^;p) ' ^^oo / MM,=o(ai,...,a^;p) 

where the second equality follows from the g-Binomial theorem [6]. 

It now suffices to show that Gg(C) = det(l + K) as in the statement of the proposition. 
From now on, all contour integrals are along Observe that we can rewrite the summation 
in the definition of /i^ so that 

C*^ 1 (mi +m2 + ■■■)! f f an A 

K^- f^miT^,... ("^i+"^2 + ■■■)! milms!--- J J f=i 

imi=k 

where w = {wi, . . . , wl) , A = (Ai, . . . , A^) and is specified by A = I'"i2™2 . . . and where the 
integrand is 

L L 



Jl(A; w; () = det 



Wiq^^ — Wj 



ll(^-q)'^e^f{wj)f{qw,)---f{q'^~' 

*j=i j=i 



The term jg multinomial coefficient and can be removed by replacing the 

inner summation by 



j ■■■ j h{n-w-0 



^1 vi^L ^■^fc.m]^ ,m2 



dwj 
27ri^ 
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with n = (ni, 

^k,mi,m2,. 

This gives 



jUl) and where 
{rii, . . . > 1 : ^ 



Hi 



k and for each j > 1, rrij of the Ui equal j}. 



L! 

L>0 ni,...,nj^>l 



ZTTl 



Now we may sum over k which removes the requirement that iii = k. This yields that 
the left-hand side of equation (14. 7p can be expressed as 

L 



L>0 ni,...,n£>l 



/.../det 


1 







J i,j=l 7 = 1 



Til — 1„ 



Wo 



) — ^. 

27ri 
(7.1) 

This is the definition of the Fredholm determinant expansion det(H-K), as desired. As these 
were purely formal manipulations, at this point to complete the proof we must justify the 
rearrangements in the above argument. In order to do this, we will show that the double sum- 
mation of (17. ip is absolutely convergent. This is the point at which the unboundedness of the 
Ca,ip contour introduces a slight divergence from the analogous proof of [16] Proposition 3.2.8 
where the contour was bounded and of finite length. 

Basically, the absolute convergence follows from the exponential decay of the function / 
as the real part of w increases to positive infinity, combined with Hadamard's inequality. Let 
us bound the absolute value of the integrand in (17. ip . Note that by assumption q^^Wi/wj — 1 
is bounded from uniformly as Wj, Wj, and rij vary. Thus, it follows that for some finite 
constant 

Since the function f{w) is bounded as w varies and has exponential decay with respect to 
the real part of u;, we can replace 

\f{w,)f{qw,) ■ ■ ■ f{q''^-'w,)\ < {B,pe-^''<^^^ 

for constants c > and B2 < oo. Thus we find that 

L 





1 


L 


det 


q^'Wi — Wj 









L>0 ■ \rt>l 



\d'w\ 



-c Rg(ui) 



(7.2) 



Since w is being integrated along Ca,^p1 the integral over w is bounded by some constant 
i?3 < 00. Finally, for \C, \ small enough the geometric series converges and it is bounded by a 
constant B^. Therefore 

^^s^{B,B,B,fL^I' 

L>0 



L! 



< 00. 



Thus we have shown that the double summation in (17. ip is absolutely convergent, completing 
the proof of Proposition 14.121 
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7.2. Proof of Theorem 14.131 This theorem and its proof are adapted from [16] Theo- 
rem 3.2.11. However, in that theorem, the u;-contour Ca.tp, was of finite length and the 
s-contour was just a vertical line. The need for slightly more involved contours comes 
from the unboundedness of the w-contour and the necessity that K(^{w,w') goes to zero 
sufficiently fast as grows along the w-contour. 

The starting point for this proof is Proposition 14.121 There are, however, two issues we 
must deal with. First, the operator in the proposition acts on a different space; second, 
the equality is only proved for |^| < for some constant C > 1. We split the proof into 
three steps. Step 1: We present a general lemma which provides an integral representation 
for an infinite sum. Step 2: Assuming C £ {C • ICI < X 4- ^^+1 derive equation (14. 9p . 
Step 3: A direct inspection of the left-hand side of that equation shows that for all C, 7^ q~^^ 
for M > the expression is well-defined and analytic. The right-hand side expression can be 
analytically extended to all C, ^ M+ and thus by uniqueness of the analytic continuation, we 
have a valid formula on all of C \ IR+ . 

Step 1: The purpose of the next lemma is to change that space we are considering and to 
replace the summation in Proposition 14.121 by a contour integral. 

Lemma 7.1. For all functions g which satisfy the conditions below, we have the identity that 

/orCG{C:|Cl<l,C^K+} 

00 

where the infinite contour Ci^2,... is a negatively oriented contour which encloses 1,2, .. . and 
no singularities of g{q^), and is defined with respect to a branch cut along z G M_. For the 
above equality to be valid the left-hand-side must converge, and the right-hand-side integral 
must be able to be approximated by integrals over a sequence of contours which enclose 
the singularities at 1,2, ... ,k and which partly coincide with Ci^2,... in- such a way that the 
integral along the symmetric difference of the contours Ci^2,... o-nd Ck goes to zero as k goes 
to infinity. 

Proof. The identity follows from Resr(— s)r(l + s) = (—1)'"'"''^. □ 

s=k 

Step 2: For this step let us assume that C ^ {C • ICI < C~^,C ^ We may rewrite 

equation fl4.8p as 

K{ni,Wi;n2,W2) = Cfl'u,i,«,2 
where g is given in equation fl4.10p . 

Writing out the M*'' term in the Fredholm expansion we have 

In order to apply Lemma mi we need to define the sequence of contours Ck (in fact we need 
only specify the contours for k large). Let Ck be composed of the union of two parts - the 
first part is the portion of the contour which lies within the ball of radius A; + l/2 centered 
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k k + 1 











1 — ; 


h 





Figure 10. Left: The contour Ck composed of the union of two parts - the 
first part is the portion of the contour which hes within the ball of radius 
+ 1/2 centered at the origin; the second part is the arc of that ball which 
causes the union to be a closed contour which encloses {1,2, . . . , k} and no 
other integers. Right: The symmetric difference between Ck and is given 
by two parts: a semi-circle arc which we call C^^'^ and a portion of i? + iM with 
magnitude exceeding k + 1/2 which we call C^^^. 



at the origin; the second part is the arc of the boundary of that ball which causes the union 
to be a closed contour which encloses {1,2,..., k} and no other integers. The contours Ck 
are oriented positively and illustrated in the left- hand- side of Figure [TOl The infinite contour 
Ci,2,... is chosen to be oriented as in the statement of the theorem (decreasing imaginary 
part). By the definition of the contours Ca,^ and we are assured that the contours Ck do 
not contain any poles beyond those of the Gamma function r(— s). This is due to the fact 
that the contours have been chosen such that as s varies, q''w stays entirely to the left of Ca,ip 
and hence does not touch w'. 

In order to apply the above lemma we must also estimate the integral along the symmetric 
difference. Identify the part of the symmetric difference given by the circular arc as C^'^^ and 
the part given by the portion of i? + iM with magnitude exceeding k + 1/2 as C^*^^ (see the 
right-hand-side of Figure [TU]). First observe that for Wi,W2 fixed, gwi,w2{Q^) stays uniformly 
bounded as s varies along these contours. Consider next (— C)*- If ~C = ^^^'^ fo^' c ^ (~^! ^) 
and r > we have (— C)* = r'^e'*"'. Writing s = x + iy we have |(— C)*| = r^e~^°". Note 
that our assumption on ( corresponds to r < 1 and a G (— vr, vr). Concerning the product of 
Gamma functions, recall Euler's Gamma reflection formula 



T(-s)T(l + s) = -^ 



Sm —TTS 



One readily confirms that for all s: dist(s, Z) > c for some c > fixed. 



TT 



sm — TTS 



c' 

< 



oTT Im(s) 
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for a fixed constant c' > which depends on c. Therefore, along the C^*^^ contour where 
s = R + iy, 

|(-C)T(-s)r(l + s)\ ~ r^e-3'"e-"l^l, 

and since a G (— tt, tt) is fixed, this product decays exponentially in \y\ and the integral goes 
to zero as k goes to infinity. Along the C^^'^ contour, the product of Gamma functions still 
behaves like c'e~'^'^' for some fixed c' > 0. Thus along this contour (again using the notation 
s = X + iy) 

|(-C)T(-s)r(l + s)\ ~ e-^V"e-"l^l. 

Since r < 1 and — (vr + cr) < these terms behave like e"'^"^^^'^'-' (c" > fixed) along the 
circular arc. Clearly, as k goes to infinity, the integrand decays exponentially in k (versus 
the linear growth of the length of the contour) and the conditions of the lemma are met. 
Applying Lemma im we find that the M^^ term in the Fredholm expansion can be written 

as 

Therefore the determinant can be written as det(l + K(^) as desired. 

Step 3: At this point we now make critical use of the choice for the contour Cq, on which 
w varies, since proving analyticity in ( of the Fredholm determinant requires the decay 
properties of the kernel with respect to w varying along Cq,^,^. 

In order to analytically extend our formula we must prove two facts. First, that the left- 
hand side of equation ( 14. 9 p is analytic for all ( ^ M+; and second, that the right-hand side 
determinant is defined (i.e., its expansion is convergent) and analytic for all ( ^ M+. 

Expand the left-hand side of equation (14.91) as 



E 

n=0 



P(A^ = n) 



(l-Cg")(l-Cg" 



where P = PMMt=o{ai,...,ajv;p)- 

Observe that for any ( ^ {'?~^^}a/=o,i,...) within a neighborhood of ( the infinite products 
are uniformly convergent and bounded away from zero. As a result the series is uniformly 
convergent in a neighborhood of any such ( which implies that its limit is analytic, as desired. 

Turning to the Fredholm determinant, we must show that the series denoted by det(l-|-i^'^) 
is an analytic function of ( away from M+. For this we will appeal to the fact that limits of 
uniformly absolutely convergent series of analytic functions are themselves analytic. Recall 
that 

It is clear from the definition of that det(A'^(wj, Wj))^^^^ is analytic in ^ away from M_|.. 
Thus any partial sum of the above series is analytic in the same domain. What remains is to 
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show that the series is uniformly absolutely convergent on any fixed neighborhood of ( not 
including R+. Towards this end consider the n^^ term in the Fredholm expansion: 




dWr. 

2m l-Ti 



dsi 
27ri 



dSr, 

27ri 



det 



N 



1 



s,)T{l + s,)i-Cr^ exp {^w,{q^^ - 1)) J] 



■m=l 



{q'^Wj/am;q)c 

[Wj/dm] q)oo 



(7.3) 



We wish to bound the absolute value of this. We may pull the absolute values inside the 
integration. Now observe the following bounds which hold uniformly over all Wj G Ca,ip and 
then all Sj G l^wj- For the first bound note that for all z : \z\ < 1, there exists a constant 
Cg such that |(2;;q')oo| < Cg. From that it follows that for any \z\ > 1, |(z;g)oo| < Cq|(z;g)fc| 
where k is such that \zq''\ < 1. This k is approximately — ln(|2r|)/ln(g) and hence bounded 
by /c < Cgln|z| for some other constant c'^. Finally |(z;g)fc| < c'^lz]'' < c'^\z\'^''^^^^^. From this 

and the fact that < 1 (recall that Re(sj) > along Vw.), it follows that for \wj\ > 1 
along Ca in we can bound 



N 

n 

m=l 



{q'JWj/am]q)c 



{wj/am] q)c 



< Ci\Wj 



lNc'^lni\w,\/d) 



(7.4) 



for some constant Ci and d = minj{aj}. For \wj\ < 1 along C^^ip, we can bound the above 
left-hand side by a constant and since \wj\ is bounded from below along Ca,ip, it follows that 
the above bound (17. 4p holds for all w G Ca^i^, possibly with a modified value of Ci. 

By Hadamard's inequality and the conditions we have imposed on we get the crude 
bound 

1 



det 



(7.5) 



.Q^^i - y i,j=i 
for some fixed constant C2 > 0. 

Finally note that by the conditions we imposed in choosing the contours, for Wj on Ca,ip 
and Sj on V^^, we have Re {wj{q^ — 1)) < where > is a constant depending on 

if G (0,7r/4). From this it follows that 

|exp {'ywj^q"^ - 1)) I < exp ( - 'yc^\wj\). (7.6) 

Taking the absolute value of (17. 3p and bringing the absolute value all the way inside the 
integrals, we find that after plugging in the results of (17. 4p . (17. 5p and (I7.6p r 



\FniC)\ < 



(ciC2)"n"/2 




-sni + s)i-cy 



\w 



Nc'ln{\w\/d) 



exp( 



-ic^\w\) 
(7.7) 



We integrate this in the s variables first. For ( ^ IR+ we would like to bound 

R we write \dz\f{z) for the integral of / along C 



®For a complex contour C and a function / : C 
with respect to the arc length \dz\. 
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/ ^|r(-.)r(i + .)(-c)^ 



for a neighborhood of C, which does not touch M+. We divide the contour of integral into two 
regions and bound the integral along each region: (1) The portion of the contour from R — id 
to 1/2 — id and then vertical to 1/2 and its reflection through the real axis; (2) The portion 
of the contour which is infinite from i? — ioo io R — id and then from R + idio R + ioo. Recall 
that by Remark 14.91 we may assume that up to constants d ~ \w\~^ and R ^ \n.\w\ for \w\ 
large enough. 

Case (1): By standard bounds |r(— s)r(l + s)\ < 1/d \w\ (since 1/ sin(x) ~ 1/x near 
X = 0). Calling r the maximal modulus over the neighborhood of \C, \ in question, it follows 
that since the s integral in Case (1) has length less than 2R (when d < 1/2), the first part of 
the integral is bounded by a constant times \w \ In |i(;|r'^3in|«>| ^j^j^ constant C3 > 0. 

Case (2): The product of Gamma functions decays exponentially in s and so the integral 
is estimated by which, by Remark 14. 9 1 is like r'^^^'^l"''. 

Summing up the above two cases we have that for \w\ large, 

/ yM|r(-s)r(l + s)(-C)1 <C4r^^'°l"lk|lnk|. (7.8) 

This estimate can be plugged in to the right-hand side of (17. 7p to reduce the bound to just 
an integral in the Wj. This integral factors and thus we have 



F„(C)|< ^^^^^V L ^kr'=^'"^'"'/'W^"^'"'k|ln|u'|exp(-7cJ^/;|) 




The integral in \w\ is clearly convergent due to the exponential decay (which easily overwhelms 
the growth of 

|y^|Arc;inH as well as the other terms). Thus the right-hand side above is 
bounded by d^n^^"^ /n\ for some constant C5. Thus -Fn(C) is absolutely convergent, uniformly 
over any fixed neighborhood of a C ^ 1^+- This implies that det(l + Kq) is analytic in C ^ R+ 
and hence completes the proof of Step 3 and hence the proof of the theorem. 

7.3. Proof of Proposition 14.201 By virtue of Lemma [8.41 it suffices to show that for some 
c,C > 0, 

\Ku{v,v')\<Ce-''\^\. (7.9) 

Before proving this let us recall from Definition 14. 81 the contours with which we are dealing. 
The variable v lies on C^^^ and hence can be written as v = a — K,cos{(p) ± ifi;sin(<y9), for 
K e M+, where the ± represents the two rays of the contour. The s variables lies on which 
depends on v and has two parts: The portion (which we have denoted I^^,c) with real part 
bounded between 1/2 and R and imaginary part ±d for d sufficiently small, and the vertical 
portion (which we have denoted T>^^\) with real part R. The condition on R implies that 
R> — Re(f ) + a + 1 and for our purposes we will assume R = — Re(f ) -|- a + 1. 

Let us denote by h{s) the integrand in (14. 2p . through which Ku{v,v') is defined. We split 
the proof into three steps. Step 1: We show the integral of h{s) over s G "Dy^^ is bounded 
by an expression with exponential decay in |f |. Step 2: We show the integral of h{s) over 
s G "D^j is bounded by an expression with exponential decay in \v\. Step 3: We show that the 
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integral of h{s) over tlie entire contour s G Vi, is bounded by a fixed constant, independent of 
V or v'. Tlie combination of these three steps imply the inequality (17. 9p and hence complete 
the proof. 

Step 1: We the various terms in h{s) separately and develop bounds for each. Let us write 
s = X + iy and note that along the contour V^^^, y G [—d, d] for d small, and x G [1/2, R]. 

Let us start with e''^'+^''/^. The norm of the above expression is bounded by the expo- 
nential of the real part of the exponent. For s along Pt,,^ 

Re(ws + = xRe{v) + — ylm{v) — —. 

We take R = — Re{v) + a + 1, d sufficiently small, and the bound Re(f ) < c' — c'\v\ for some 
constants c',c' (depending on (p), to deduce 

Re(f s + < c — c\v\x 

for some constants c, c > 0, from which 

Let us now turn to the other terms in h{s). We have 

and we may also bound 
T(v — a 



"m J 



T(s + V — ttr. 





1 


<c, 




V + s — v' 



<a ir(-s)r(i + s)i < c, 



for some constants C > (which may be different in each case). The first bound comes from 
the functional equation for the Gamma function, and the last from the fact that s is bounded 
away from Z. 

Combining these together shows that for large, the portion of the integral of h{s) for s 
in □ is bounded by (recall s = x + iy) 

^-CT\v\x+xln\u\-arg{u)y ^ Ce~'^^^^ 



for some constants c,C > 0. Since for \y\ in a bounded set, everything starting from the 
integration path is clearly bounded, the bound holds. 

Step 2: As above, we consider the various terms in h{s) separately and develop bounds for 
each. Let us write s = R + iy and note that s G corresponds to y varying over all |y| > ci. 
As in Step 1, the most important bound will be that of e''^*+'^'*^/2. 
Observe that 

Observe that because if> G (0, 7r/4) and R = — Re(f ) + a + 1, 

—^ + — + Re{v)R<c-c\v\^ (7.10) 
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for some constants c, c > 0. Thus 



Re(.. + sV2) < -1^±M!^ + c-c\v\\ (7.11) 



Let us now turn to the other terms in h{s). We bound 



s\ ^ ^R\n\u\-y a.Tg{u) 



By standard bounds for the large imaginary part behavior we can show 



r(s + t> — a„ 

for some constant C > sufficiently large. Also, \l/{y + s — v')\ < C for a fixed constant. 
Finally, the term 

|r(-s)r(i + s)| < Ce-"i^i, 

for some constant C > 0. 

Combining these together shows that the integral of h{s) over s in is bounded by a 
constant time 

/"exp /^_^(^JlI^^M]_ _ Tc\v\^ + R\n\u\ - y arg{u) - n\y\ + N^\y\) dy. (7.12) 

We can factor out the terms above which do not depend on y, giving 

exp [^—Tc\v\ + Kln\u\) / exp I — r y a.Tg[u) + JM — \y\ \ dy. 

V 2 2 / 

Notice that the prefactors on y and |y| in the integrand's exponential are fixed constants. We 
can therefore use the following bound that for a fixed and 6 G M, there exists a constant C 
such that 

f e-^^y+^^'+^\y\dy < Cel'^^l, (3 > 0. 
Jr 

Using this we find that we can upper-bound fl7.12p by 

exp (-rc|w|^ + R\n \u\ + c'\v\) . (7.13) 

For \v\ large enough the Gaussian decay in the above bound dominates, and hence integral 
of h{s) over s in | is bounded by 

fj^-c\v\ 

for some constants c,C > 0. 

Step 3: Since v' only comes in to the term l/{v + s — v') in the integrand, it is clear that the 
above arguments imply that the integral of h{s) over the entire contour s eV^j is bounded by 
a fixed constant, independent of v or v'. This completes the third step and hence completes 
the proof of Proposition 14.201 
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7.4. Proof of Proposition I47l9l Fix ?7,r > 0. We are presently considering tlie Fredliolm 
determinant of the kernels K^^ and Ku restricted to the fixed finite contour Ca,ip,<r- By 
Lemma 18.51 we need only show convergence as e — )■ of the kernel Kl^{v,v') to Kuiv,v'), 
uniformly in v,v' G Ca,ip,<r- This is achieved via the following lemma. 

Lemma 7.2. For all rj' > there exists Eq > such that for all e < Eq and for all 

\K:^{v,v')-K^{v,v')\<ri'. (7.14) 

Proof. The kernels and Ku are both defined via integrals over s. The contour on which s 
is integrated can be fixed for {e < Eq) to equal Vy, which is the s contour used to define Ku- 
The fact that the s contours are the same for K^ and Ku is convenient. The proof of this 
lemma will follow from three claims. The first deals with the uniformity of convergence of 
the integrand defining K^ to the integrand defining Ku for s restricted to any fixed compact 
set. 

Before stating this claim, let us define some notation. 

Definition 7.3. Let T>y y^j = {s G P^, : \s\ > M} be the portion ofV^ of magnitude greater 
than M and similarly let Vy^^M = {s & Vy : \s\ < M}. Let us assume M is large enough 
so that T>y^yM is the union of two vertical rays with fixed real part R = — Re(w) + a + 1. 
Assuming this, we will write s = R + iy. Then for yM = {M"^ — (— Re(w) + a + 1)^)""^/^, the 
contour Vy^^M = {R + iy : \y\ > yAi}- 

Claim 7.4. For all rj" > and M > there exists Eq > such that for all e < Eq, for all 

V, v' G Ca,^,<r, o,nd for all s G Vy^^M, 



h''{s)-T{-s)T{l + s)ll ^^"^ 



, r(s + V - am) V + S 

m=l ^ ' 



< (7.15) 



where h'^ is given in Ili4-i3{ ). 



Proof. This is a strengthened version of the pointwise convergence in ( I4.14p through fl4.17p . It 
follows from the uniform convergence of the Tq function to the F function on compact regions 
away from the poles, as well as standard Taylor series estimates. The choice of contours is 
such that the pole arising from l/[v + s — v') is uniformly avoided in the limiting procedure 
as well. □ 

It remains to show that for M large enough, the integrals defining K^{v, v') and Ku{v, v') 
restricted to s in Vy^-^M, have negligible contribution to the kernel, uniformly over v,v' and 
E. This must be done separately for each of the kernels and hence requires two claims. 

Claim 7.5. For all rj' > there exists Mq > and Eo > such that for all e < Eq, for all 

V, v' G Ca,^^<^r, and for all M > Mq, 

I dsh\s) <ri'. 



FREE ENERGY FLUCTUATIONS FOR DIRECTED POLYMERS 



63 



Proof. We will use the notation introduced in Definition 17.31 and assume Mq is large enough 
so that Vy^yM is only comprised of two vertical rays. 

Let us first consider the behavior of the left-hand side of fl4.17p . The magnitude of this 
term is bounded by the exponential of 

Re{Te~'s + e~^Tq"{q' - 1)). 

This equation is periodic in y (recall s = R+iy) with a fundamental domain y G [—7re~^, 7re~^]. 
For e~^7r > \y\ > y^ for some y^ which can be chosen uniformly in v and e, the following 
inequality holds 

Re(TC-^s + - 1)) < -Ty'^/% 

This can is proved by careful Taylor series estimation and the inequality that for x G [— tt, tt], 
cos(x) — 1 < — This provides Gaussian decay in the fundamental domain of y. 

Turning to the ratio of g-Gamma functions in fl4.16l) . observe that away from its poles, the 
denominator 

1 ,c'rw 



Kce'^^') (7.16) 



where c, c' are positive constants independent of e and v (as it varies in its compact contour) 
and /^(s) = dist(Im(s), 27r£~^Z). This establishes a periodic bound on this denominator, 
which grows at most exponentially in the fundamental domain. The numerator Tq{v — am) 
in fl4.16p is bounded uniformly by a constant. This is because the v contour was chosen to 
avoid the poles of the Gamma function, and the convergence of the g-Gamma function to the 
Gamma function is uniform on compact sets away from those poles. 

Finally, the magnitude of (14.141) corresponds to \u^\ and behaves like e-^in|«|+j/arg(M)_ ^hus, 
we have established the following inequality which is uniform in v, v' and e as, y varies: 



< c'e'-^'^'^'/^+'^'^l-^'^'^l (7.17) 



il-q)^) q'+-~q-' A i r,(s + - lng(5^)) 

for some constant c" > 0. Notice that this inequality is periodic with respect to the funda- 
mental domain for y G [—7rs^^,7rs~^]. 

The last term to consider is r(— s)r(l + s) which is not periodic in y and decays like e~'^'^' 
for ?/ G M. Since 'P„,>m is only comprised of two vertical rays we must control the integral of 
h'^{s) for s = R + iy and \y\ > yu- By making sure M is large enough, we can use the periodic 
bound (17.171) to show that the integral over yu < \y\ < e'^ii is less than 77 (with the desired 
uniformity in and e. For the integral over \y\ > e^^ir, we can use the above exponential 
decay of r(— s)r(l + s). On shifts by 27re~^Z of the fundamental domain, the exponential 
decay of r(— s)r(l + s) can be compared to the boundedness of the other terms (which is 
certainly true considering the bounds we established above). The integral of each shift can 
be bounded by a term in a convergent geometric series. Taking Eq small then implies that 
the sum can be bounded by rj' as well. Since 77' was arbitrary the proof is complete. □ 

Claim 7.6. For all t]' > there exists Mq > such that for all v, v' G Ca^^p^^T, and for all 
M > Mo, 

^ ^(v - a™) M^e^'"^+"'/2 



f dsT{-s)T{l + s)Y[ 



< rj 
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Proof. The desired decay here comes easily from the behavior of vs + as s varies along 
Vjj^yM- As before, assume that Mq is large enough so that this contour is only comprised of 
two vertical rays and set s = i?+i?/ for ?/ G M for \y\ > Um- As in the proof of Proposition 14.201 
given in Section 17.31 one shows that 

uniformly over w, v' G Ca^^^<R, and for all M > Mq. This behavior should be compared to that 



r{v-am) 



of the other terms: |r(-s)r(l + s)| ^ e"^!*'!: \u'\ = g-^i^l^l+f'^-'gH; ^, ^ , 

I \ ^ \ ' ^1 111 ) r(s+v-amj 

and 1 1/ + s + ?;') I < C as well. Combining these observations we see that the integral decays 
in \y\ at worst like Ce~'^y^~^'^'^^^ . Thus, by choosing M large enough so that um 3> 1 we can be 
assured that the integral over \y\ > Um is as small as desired, proving the above claim. □ 

Let us now combine the above three claims to finish the proof of the Proposition 14.191 
Choose 7]' = ri/3 and fix Mq and as specified by the second and third of the above claims. 
Fix some M > Mq and let L equal the length of the finite contour Vv,<m- Set t]" = ^ and 
apply Claim [731 This yields an Eq (which we can assume is less than ^q) so that f l7.15p holds. 
This implies that for e < Eq, and for all v, v' G Co,,^,<.r, 

From the triangle inequality and the three factors of rj'/S we arrive at the claimed result of 
f l7.14p and thus complete the proof of the lemma and hence also Proposition 14.191 □ 



7.5. Proof of Proposition I47l8l The proof of this proposition is essentially a finite e (recall 
q = e~^) perturbation of the proof of Proposition 14.201 given in Section 17.31 The estimates 
presently are a little more involved since the functions involved are g-deformations of classic 
functions. However, by careful Taylor approximation with remainder estimates, all estimates 
can be carefully shown. By virtue of Lemma [8.41 it suffices to show that for some c, C > 0, 

\K{v,v')\<Ce-'\''\. (7.18) 

Before proving this let us recall from Definition l4.16l the contours with which we are dealing. 
The variable v lies on C^^^ V ^ (0,vr/4). The s variables lies on P„ which depends on v 

and has two parts: The portion (which we have denoted 2^i,,iz) with real part bounded between 
1/2 and R and imaginary part ±d for d sufficiently small, and the vertical portion (which we 
have denoted I^i,,iz) with real part R. The condition on R is that R> — Re(f ) + a + 1 and 
for our purposes we can take that to be an equality. 

Let us recall that the integrand in fl4.12p . through which K^{v,v') is defined, is denoted 
by h'^{s). We split the proof into three steps. Step 1: We show the integral of /i'^(s) over 
s G T>^^^ is bounded for all e < Sohy an expression with exponential decay in Step 2: 
We show the integral of h'^{s) over s G is bounded for all e < e^hy an expression with 
exponential decay in \v\. Step 3: We show that for all e < Eq, the integral of h'^{s) over the 
entire contour P„ is bounded by a fixed constant, independent of v or v'. The combination 
of these three steps imply the inequality fl7.18p and hence complete the proof. 
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Step 1: We consider the various terms in h'^{s) separately (in particular we consider the left- 
hand sides of (14.141) through (14.171) ) and develop bounds for each which are valid uniformly 
for e < Eq and e < small enough. Let us write s = x + iy and note that along the contour 
1^v,\zy y ^ [~d, d\ for d small, and x G [1/2, R]. 

Let us start with the left-hand side of (14.171) which can be rewritten as 

exp (r Re{e-^s + e~^q\q' - 1))) . 

The norm of the above expression is bounded by the exponential of the real part of the 
exponent. For G (0, 7r/4) one shows (as a perturbation of the analogous estimate in Step 1 
of the Proof of Proposition 14.201) via Taylor expansion with remainder estimates that 

r Re(£:~^s + e~^g''(g'' - 1)) < c - tc\v\x, 

for some constants c, c. The above bound implies 

|exp (r Re(e-is + e-^q'^iq' - 1))) | < Ce'^l'^l^ 

Let us now turn to the other terms in h'^{s). We bound the left-hand side of (I4.14p as 

-c 



N 



We may also bound the left-hand sides of (I4.15p and (I4.16p . as well as the remaining product 
of Gamma functions by constants: 



Tqi^v - Ing(am)) 



q" In q 



q- 



ir(-s)r(i 



for some constants C > (which may be different in each case). The first bound comes 
from the functional equation for the g-Gamma function, and the last from the fact that s is 
bounded away from Z. 

Combining these together shows that for \v\ large. 



[ h''{s)ds 


iv 







for some constants c', C > 0, while for bounded \v\ the integral is just bounded as well. 

Step 2: As above, we consider the various terms in h'^{s) separately and develop bounds for 
each. Let us write s = R + iy and note that s G corresponds to y varying over all 
\y\ > d. Three of the terms we consider (corresponding to the left-hand sides of (I4.15p . (I4.16P 
and (I4.17P ) are periodic functions in y with fundamental domain y G [—ne~^, ne^^. We will 
first develop bounds on these three terms in this fundamental domain, and then turn to the 
non-periodic terms. 

We start by controlling the behavior of the left-hand side of (I4.17P as y varies in its 
fundamental domain. For each (f < 7r/4 there exists a sufficiently small (yet positive) constant 
c' such that as y varies in its fundamental domain 

rRe(e-^s + £-V(g' - 1)) < c'rRe{vs + /2). 
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On account of this, we can use the bound ( 17. lip from the proof of Proposition 14.201 This 
imphes that 

{y + lm{v)Y 



rRe(e"^s + £"V(g' - 1)) < cV 



— c\v\ 



Let us now turn to the other y-periodic terms in h'^{s). By bounds for the large imaginary 
part behavior of the g-Gamma function we can show 



Tg{v - Ing(am)) 



for some constants c, C > where /^(s) = dist(Im(s), 27r£:"-'^Z). Note that as opposed to 
f l7.16p when \v\ was bounded, in the above inequality we write f^{s + v) in the exponential 
on the right-hand side. This is because we are presently considering unbounded ranges for v. 
Also, we can bound 

In q 



< c 



for some constant C > 0. 

The parts of h'^{s) which are not periodic in y can easily be bounded. We bound the 
left-hand side of f l4.14p as in Step 1 by 



-c 



Finally, the term 



{i-qr 

|r(-s)r(i + < Ce-^\y\ 



for some constant C > 0. 

We may now combine the estimates above. The idea is to first prove that the integral on the 
fundamental domain y G [— vre"^, vre"^] is exponentially small in \v\. Then, by using the decay 
of the two non-periodic terms above, we can get a similar bound for the integral as y varies over 
all of M. For j G Z, define the j shifted fundamental domain as Dj = je~^2'K + [—e~^Ti^ e'^rc]. 
Let 



/, := / h\R + iy)dy 



and observe that combining all of the bounds developed above, we have that 



|/,| < C 



where 



Fi{y) = exp cV 



{y + lm{v)Y 



F,{y)F,{y)dy, 



c\v\ \ + c" f'^{s + v) + x\n\u\ 



F2{y) = exp (-{y + je ^2%) aTg{u) - TT\y + je ^2-r\) . 

The term Fi{y) is from the periodic bounds while -^2(2/) from the non-periodic terms (hence 
explaining the je~^27r shift in y). By assumption on u, we have — arg(n) — tt = 5 < c for 
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some 6. Therefore 

form some constants c,C > 0. Thus 

Just as in the end of Step 2 in the proof of Proposition I4.2UI we can estimate the integral 



for some constants C, c > 0. This imphes 
Finally, observe that 



h''{s)ds 



where C is independent of e as long a.s e < Eq for some fixed Sq- This is the bound desired 
to complete this step. 

Step 3: Since v' only comes in to the term ^s+vl'^v' the integrand, it is clear that the above 

arguments imply that the integral of h'^{s) over the entire contour s G I5„ is bounded by a 
fixed constant, independent of v or v'. This completes the third step and hence completes 
the proof of Proposition 14.181 

8. Appendix 

8.1. Two probability lemmas. 

Lemma 8.1 (Lemma 4.1.38 of [IS]). Consider a sequence of functions {/n}n>i niapping 
M — [0,1] such that for each n, fn{x) is strictly decreasing in x with a limit ofl at x = —oo 
and at X = oo, and for each 6 > 0, on ]R\ [—6, 6], fn converges uniformly to lx<o- Consider 
a sequence of random variables Xn such that for each r G M, 

E[/„(X„-r)]^p(r) 

and assume that p{r) is a continuous probability distribution function. Then Xn con- 
verges weakly in distribution to a random variable X which is distributed according to 
P(X < r) =p{r). 

Lemma 8.2 (Lemma 4.1.39 of [IS]). Consider a sequence of functions {/n}n>i mapping 
M — 7- [0, 1] such that for each n, fn{x) is strictly decreasing in x with a limit ofl at x = — oo 
and at X = oo, and fn converges uniformly on M. to f . Consider a sequence of random 
variables Xn converging weakly in distribution to X . Then 

E[fn{Xn)] ^ E[/(X)]. 
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8.2. Some properties of Fredholm determinants. We give some important properties 
for Fredholm determinants. For a more complete treatment of this theory see, for exam- 
ple, [59]. 

Lemma 8.3 (Proposition 1 of ^S])- Suppose t Tf is a deformation of closed curves and a 
kernel L{ri,ri') is analytic in a neighborhood ofTf x C for each t. Then the Fredholm 
determinant of L acting on Ft is independent oft. 

Lemma 8.4. Consider the Fredholm determinant det(l + i^)i2(r) on an infinite complex 
contour T and an integral operator K onV. Parameterize T by arc length with some fixed 
point corresponding to r(0). Assume that \K{v,v')\ < C for some constant C and for all 
v,v' eT and that either of the following two exponential decay conditions holds: There exists 
constants c, C > such that 

\K{T{s)J{s'))\<Ce-'^\^\ 

Then the Fredholm series defining det(l + _R')i2(r) is well-defined. Moreover, for any r] > 
there exists an > such that for all r > Tq 

I det(l + K)L2(r) - det(l + /r)L2{r.)l < V 
where Tr = {r(s) : |s| < r}. 

Proof. The Fredholm series expansion (14. 6 p is given by 

det(l + i^)z.2(r) = [ dsi- - [ dsndet{K{T{si),T{sjm^^, (8.1) 

n>o -^r 

is well-defined since by using Hadamard's bouncil one gets that 

n 

\det{K{r{s,),r{sj)))lj^^\ < JJe-'^l^'^l (8.2) 

i=i 

which is absolutely integrable / summable. To show is det(l + K)l2(y^-^ — det(l + K)L2(^r) 
as r — )■ oo. From (18. ip one immediately gets that 

det(l + K)L^r.) = det(l + PrK)L2(r). 

where is the projection onto F^. The kernel {PrK){si, sj) converges pointwise to K{si, Sj) 
and (18. 2 p provides a in r uniform, integrable / summable bound for det{K{T{si), T{sj)))^j^i. 
Therefore, by dominated convergence as r — > oo the two Fredholm determinant converge. □ 

Lemma 8.5. Consider a finite length complex contour T and a sequence of integral operators 
on r, as well as an addition integral operator K also on T . Assume that for all t] > 
there exists Eq such that for all e < Eq and all z, z' G T, \K'^{z, z') — K{z, z')\ < rj and that 
there is some constant C such that \K{z, z')\ < C for all z, z! G F. Then 

limdet(l + i^^)L2{r) = det(l + K)i^nj^). 



Hadamard's bound: the determinant of a n x ti matrix with entries of absolute value not exceeding 1 is 
bounded by n"/^. 
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Proof. As in Lemma [83] one writes the Fredholm series. Since F is finite, the Fredholm deter- 
minants det(l + K)i2(j-j is well-defined because \K{z, z')\ < C (use Hadamard's bound). By 
assumption, K"^ converges pointwise to K and we have the uniform bound \K'^{z, z')\ < C + rj. 
This ensures that det(l -|- K^)L2(r^ is well-defined and that we can take the limit inside the 
Fredholm series, providing our result. □ 

8.3. Reformulation of Fredholm determinants. 

Lemma 8.6. Let Kai in / 1 5. 5\) . as in 1^5. 4\ ), and Kai the Airy kernel. Then it holds 

det(l + ^Ai)L2(C») = det(l - i^Ai)L2(r,oo)- 

Proof. The integration path in (15. 5p can be chosen to have R.e{z) > and since Re(w;) < 
for w E Cw, 'we can use 



z — w 



to write KAi{w,w') = -{AB){w,w') with A : L'^{C^) L'^0^+) and B : L^{R+) L'^iC^) 
have kernels 

A{W,\) = e-V3+«.(r+A)^ ^ / 

We also have 

{BA){r],r]') = ^ / dwB{T],w)A{w,r]') 

= -, TTT / dw \ dz ^ — , , = KpArt ^ r,r\ ^ r\. 

Then, since det(l — AB')i2(q^'s^ = det(l — BA)i2(jg_^^ = det(l — -ft'Ai)L2{r,oo) we get the claimed 
result. The first equality is a general result which applies as long as AB and BA are both 
trace-class operators [59]. □ 

Lemma 8.7. Let Kbbp as in / 1 5. 23\) . as in Theorem \2.1\ (b), and Kbbp as in (T^. Then 
it holds 

det(l + Kbbp)l2{c^) = det(l - KBBp)LHr,oo)- 

Proof. The proof is as the one of Lemma [8T6| except that in A{w, A) is multiplied by YYk=i lijhr 
and 5(A, w') by UZii^ -h). " □ 



Lemma 8.8. Let -R'cdrp as in ^6. 61) . as in ( fg. and A'cdrp the CDRP kernel given in 
( [i. 7p . Then it holds 

det(l + -K^cdrp)l2(c„) = det(l - -ft'cDRp)L2(iR^). 

Proof. Using 
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we get 



KcDRp{w,w') = / dr]A{w,r])B{r],w') 
Jr+ 

with B{ri,w') = e^'^' and 



A{w,r]) = — / dz—— — - ^\ 

2m . sm{7i{z - w]a] 



4(T 



Thus det(l + -ft'cDRp)L2(c„) = det(l - Kcdrp)l2{r+) where Kcdrp = -BA, namely 

i^CDRp(?7, V) = -TT-T / dwB{r],w)A{w,r]') 
Zm y_ 1 +iR 



aw / a;2- 



(8.3) 



(27ri)2 y_ 1 7 1 sin(7r(2; - e-'/3— ' 

The next step uses the following identity: for < Ke{u) < 1 it holds 

dt 



dt. 



sin(7rM) Jjg S" + e* 
from which, for < Ke{u) < l/a it holds 



sml^vrcTMj 

We can use wit u = z — w and obtain 



S + e-*/- 



i^CDRp(r/,V) = / dt—^(^ [ dwe-'/'^-^^^A [ dze^'^'-^^^'^A 

= / dt ^J_,, Ai{ri + t)Ai{ri' + t), 

Jr o -\- e ' 

(8.4) 

the expression of (11.71) . □ 
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